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Abstract 

We present a new energy-dependent partial-wave analysis of all antiproton-proton elastic {pp — )• pp) 
and charge-exchange {pp — )• nn) scattering data below 925 MeV/c antiproton laboratory momen- 
tum. The long-range parts of the chiral one- and two-pion exchange interactions are included 
I exactly. The short-range interactions, including the coupling to the mesonic annihilation chan- 

Qj" nels, are parametrized by a complex boundary condition at a radius of r = 1.2 fm. The updated 

■ database, which includes significantly more high-quality charge-exchange data, contains 3749 scat- 

tering data. The fit results in Xmm/^df = 1-048, where A'^df = 3578 is the number of degrees of 
freedom. We discuss the description of the experimental data and we present the antiproton-proton 
phase-shift parameters. 

PACS numbers: 13.75.Cs, ll.80.Et, 12.39.Fe, 21.30.Cb 
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I. INTRODUCTION 



The antinucleon-nucleon (NN) interaction at low energies is of fundamental interest, but 
progress towards understanding it has always been hindered by the lack of scattering data. 
Major steps forward were taken at the Low Energy Antiproton Ring (LEAR) at CERN 
in the 1980 's and the early 1990 's. For the first time, good-quality data became available 
for the total cross section and the total annihilation cross section as function of antiproton 
laboratory momentum (piab), for the analyzing power in antiproton-proton elastic scattering 
[pp — pp), and for the differential cross section and analyzing power in charge-exchange 
scattering {pp — > nn), at antiproton momenta above about 200 MeV/c. Unfortunately, 
LEAR was closed in 1996 and pp scattering experiments came to a halt. However, the 
enormous physics potential of a low-energy antiproton beam is clear, especially when it can 
be polarized, and in recent years the interest to investigate pp scattering has been revived, 
for instance by the collaboration for Polarized Antiproton experiments (PAX) Ij. 

The dominant feature of antiproton-proton scattering at low energy is the annihilation 
into mesons, a complex multiparticle process that is difficult to model. In pre-LEAR days, 
some qualitative understanding was obtained by using simplified prescriptions, such as a 
simple absorptive boundary condition 2|-4| or a state-independent two- or three-parameter 
optical potential j5-12|. These models could describe the integrated total, annihilation, 
and charge-exchange cross sections, but not the differential observables. Motivated by the 
experiments at LEAR, more sophisticated NN models were developed in order to attempt 



a more quantitative fit to the data. Examples are the Paris optical-potential model |l3l4l8| 



and the Nijmegen |l9l. |20|| and Pittsburgh [21] coupled-channels models. 

In Refs. 22l-l24| an energy- dependent partial-wave analysis (PWA) of all pp scattering 
data below piab = 925 MeV/c was developed, in order to arrive at a model-independent 



description of the AW interaction. The met 



PWAs of the pp and np scattering data 



lod of analysis was adapted from the Nijmegen 
25|-|29| . These PWAs exploit as much as possible 



our knowledge about the interaction in the description of the energy dependence of the scat- 
tering amplitudes. The long-range interactions, which are responsible for the rapid energy 
variations of the amplitudes, are included exactly in the Schrodinger equation, while the slow 
energy variations due to the essentially unknown short-range interactions are parametrized 
phenomenologically by a state- and energy-dependent boundary condition at some radius 



b. In this way, an economic and model- independent high-quality description of the 



scattering database is possible. In the NN case |22N24| . one assumes that the long-range po- 
tential is given by the charge-conjugated version of a corresponding nucleon-nucleon {NN) 
potential, and, by implementing a complex boundary condition, one bypasses with this 
strategy as well our lack of knowledge of the short-range annihilation dynamics. 



There are two important reasons to update the pp PWA of Ref. 23|]. The first and 
perhaps main motivation is the renewed experimental interest in NN scattering. The second 
reason is theoretical and is motivated by the progress reached in the last two decades in the 
understanding of the AW interaction within the framework of chiral effective field theory. In 
particular, the pp and np PWAs have been updated by including, next to the electromagnetic 
and the one-pion exchange (OPE) potential, the long-range parts of the chiral two-pion 



exchange 



fT 



potential 30 



) potential 28|, |29|, instead of the heavy-boson exchanges of the Nijmegen 



31| . thereby 



improving even more the model independence and the quality 



of the AW PWAs of Refs. 



25 



27|. Motivated by that success, we include here as well the 



charge-conjugated TPE potential in the long-range WW interaction, instead of the charge- 
conjugated heavy-boson exchanges that were used in Ref. [23 1. 

At the same time, we take the opportunity to update the database of pp scattering data. 
The database constructed in Ref. [23j included all scattering data published in a regular 
physics journal up to early 1993. A number of high-quality data sets from LEAR became 
available only later, in particular differential cross sections and analyzing powers for the 
charge-exchange reaction pp — )■ nn. Also the first measurements of the depolarization and 
spin-transfer observables for pp — )■ nn were published only later. These data sets can be 
included now and they provide significant new constraints on the PWA solution. 

The organization of our paper is as follows: In Sec. [IT] the method of PWA developed 



in Ref. 



23| is reviewed. We summarize only the main points in order to make this paper 



self-contained and we emphasize the differences of our PWA with Ref. [23|. In Sec. IIIII we 
discuss the boundary condition that parametrizes the short-range interaction. In Sec. llVl the 
long-range WW potential is discussed, in particular its chiral TPE component. In Sec. |Vl the 
new database is discussed and the statistical methods are reviewed. In Sec. |VT] we present 
the results of the PWA and discuss the description of the measured observables. In Sec. IVIII 
we present the WW S matrix and phase-shift parameters. We conclude in Sec. IVIIII An 
Appendix is devoted to a study of the statistical quality of the database. 



II. THE METHOD OF ANALYSIS 



For states with total angular momentum J, the radial part of the wave function for 
the antiproton-proton system, $'^(r), is obtained by solving the coupled-channels radial 
Schrodinger equation 

$'^(r) = , (1) 

which is a differential equation in channel space. We include the channels pp and nn. It is 
important to use this physical basis instead of the isospin basis, in order to be able to include 
the long-range electromagnetic interactions and to treat the threshold for charge-exchange 
scattering pp — )• nn at piab — 99 MeV/c (or Tiab — 5.2 MeV) properly, which gives a much 
better description of the low-energy charge-exchange data. In Eq. ([1]), p is a diagonal matrix 
with the channel momenta Pa in the center-of-mass system, m is a diagonal matrix with the 
reduced mass of the two scattered particles in channel a (so nia = Mp/2 or M„/2), and 
V"^ is the potential with matrix elements {£' s'a'lV"^ {r)\i s a) . For partial waves with i = J, 
s = 0, 1, or £ = 1, J = 0, the matrices are 2x2, and for partial waves with i = J ± 1 
(<^ > 1); s = 1, coupled by the tensor force, the matrices are 4x4. The relation between 
the total energy ^/s in the center-of-mass system and the channel momentum is given by 
the relativistic expression \s = pl + '^m?^. 

We solve Eq. ([T]) numerically, starting with the boundary condition at r = 6, up to 
"r = oo," which in practice is a point outside of the range of the strong interaction. The 
asymptotic form of $"^(r) for r = oo can be written as 

KirV-^ ^[H,{pr)S' + H,{pr)\ , (2) 

where 5*"^ is the partial-wave 5* matrix and Hi and H2 are diagonal matrices. For the pp 
channel, where the Coulomb force acts, the entries are given by 

Hf\r],pr) = F(,{ri,pr) - iGi{ri,pr) , 

H} '{r],pr) = Fi{r],pr) + iGe{r],pr) , (3) 

where Fi and Gi are the standard regular and irregular Coulomb wave functions; 77 = a/wiab 
is the relativistic Coulomb parameter, where a is the fine-structure constant and fiab is the 
velocity of the incoming antiproton in the laboratory frame. The asymptotic behavior of F^ 
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and Ge is 



Fi{ri,pr) sin 
Gi{r],pr) 



pr 



- + o-^ -r]\n{2pr) 



r— >oo 
~ COS 



pr — £— + a£ — r) ln(2pr) 



(4) 



where the Coulomb phase shift is ae = arg r{i + 1 + irj). For the nn channel, rj = 0, and 



Ge{0,p) = -pneXp) 



(5) 



where ji{p) and ni{p) are the spherical Bessel and Neumann functions. The S matrix is 
obtained from the matching condition 



(6) 



where is the numerical solution of Eq. ([T]) and is given by Eq. ([2]). The Wronskian 
is defined by 

_ 1 _ 1 

(7) 



$2) = $?-$2 - ^7-^2 

m m 



where the prime denotes differentiation with respect to r and "T" means transposition. This 
gives for the partial-wave 5* matrix 



= - 



V m 



'mp 



Jmp V ^ 



(8) 



the prime on the Hankel functions denotes differentiation with respect to the argument pr. 
Due to the presence of the long-range electromagnetic interaction, care has to be taken to 



23 



25 



32| . We include in the potential 



define the S matrix [i.e. the phase-shift parameters) 
the long-range parts of the Coulomb, the magnetic- moment, and the strong (one- and two- 
pion exchange) interactions, V = Vc + Vmm + Vn- We integrate the Schrodinger equation up 
to a point outside the range of the strong interaction, where we match to Coulomb (for pp) 
and Bessel (for nn) wave functions. The S matrix is therefore defined with respect to the 
Coulomb force that acts in the pp channel. Because we need to include the infinite-range 
Coulomb interaction and part of the magnetic-moment interaction in all partial waves, but 
the finite-range nuclear interaction only up to some maximum value of J, we decompose the 
S matrix in order to split off the Coulomb part and the magnetic-moment part as 

'c+mm) K'^c+mm+n 
5 



Sc+MM+N — 1 — [Sc — 1) + Sj^^ [SS+MM ~ 1) ^c'^ + 



q1/2 fqC 



,1/2 ^ qC+MM 



1) {^c+mm) 



1/2 r-1/2 



(9) 



where Sc is the Coulomb S matrix with matrix elements {£'s'\Sc\(^ s) = Su'Sss' exp(2z(T^) 
in the pp channel and zero in the nn channel. In Eq. (Q we used matrix notation, because 
the magnetic-moment interaction contains a tensor part and the S matrix is not diagonal 
in orbital angular momentum; its square root is well-defined, however. 
The scattering amplitude is correspondingly decomposed as 

Mc+mmMO) = Mcie) + MS^maM + M^X^^^^id) , (10) 

where MciO) is the Coulomb scattering amplitude, M^j^^^^jiO) is the magnetic-moment 
scattering amplitude in the presence of the Coulomb interaction, and ^c+MM+iv(^) 
scattering amplitude for the strong interaction in the presence of the Coulomb and magnetic- 
moment interactions. The matrix elements of Mc{9) for pp scattering are given by 

p\l — cos U) 

= YpJ^Tey^v ■ (11) 
The matrix elements ^c+MAf(^) magnetic- moment interaction are calculated in 

Coulomb distorted- wave Born approximation {23!, 32 1. The partial- wave decomposition of 
the nuclear scattering amplitude is given by 



{s'm'a' I M^Xmm+n (^) k m a) = VM2£+1) C'.t^i C i y^.^, (6) 

ee'j 

{£sa\sj {Sc+mm) {^c+mm+n ^ ^) {^c+mm) \^ ^ /'^Wa , (12) 



where a denotes the channel pp or nn. Because S'^Xmm+n difficult to calculate it 



IS 



approximated by S^Xmm+n — ^c+ni where S'^+^y is the S matrix for the strong interaction 
in the presence of the Coulomb interaction. From the scattering amplitude on the spin- 
singlet, spin-triplet basis, all the observables can be calculated [33|. 



III. THE BOUNDARY CONDITION APPROACH 



The coupled-channels Schrodinger equation, Eq. ([T]), is solved with a boundary condition 
at a radius r = b, for each energy and for each partial wave. The fit to the data is not very 
sensitive to the exact value of b, but in our case an optimal value b = 1.2 fm was found. For 
the specific form of the partial- wave boundary condition we define the P matrix 34, |35| by 



($ 



dr 



(13) 



- r=b 



where $'^(r) is the radial wave function. The P matrix parametrizes the comphcated short- 
range interaction of the pp system. The couphng of the pp and nn channels to the mesonic 
annihilation channels is taken into account by a complex P matrix. 

The P matrix is a powerful tool in a PWA, since it provides the separation between the 
long-range interaction, which is relatively model independent and taken into account exactly 
in the Schrodinger equation, and the short-range interaction, which is essentially unknown 
and parametrized completely phenomenologically. The long-range interactions cause the 
rapid energy dependence of the scattering amplitudes, while the short-range interactions 
result in slow energy variations. The results, for that reason, do not depend on the details of 
the short-range interactions. We therefore choose a simple parametrization for the P matrix, 
which corresponds to a state-dependent, i.e. spin- and isospin-dependent, short-range optical 
potential. We assume that the interaction in each partial wave can be parametrized by 
a complex spherical well, the depth of which is different for elastic and charge-exchange 
scattering, i.e. for / = and 1 = 1. For a single- channel partial wave with orbital angular 
momentum i, isospin /, and with the spherical well V/ + iWj, the P matrix is given by 

Pi = p'hJ[{p'h)/Jt{ph) , (14) 

where Jt{p) = pjt{p) and p'^ = p2 _ 'M{y^ + ^p^^)^ ^here M = {Mp + M„)/2. 

The P matrix is calculated on the isospin basis and then transformed to the physical 
particle basis with the channels pp, nn. For the uncoupled partial waves with i = J, 
s = 0,1, OT £ = 1, J = 0, itis therefore a 2x2 matrix. For the partial waves with i = J ±1 
(J > 1), s = 1, coupled by the tensor force, we introduce for each value of the isospin I an 
additional mixing angle 6ij between the partial waves with i = J — 1 and i = J + 1. We 
write 

cos 6'/j sin 6^77 \ ( i ) I cos 9 u -smOu 
-sin 6*77 cos 6jj J y Pj+iJ [sin Ojj cos 6ij 

where Pj-i and Pj+i are the single-channel P matrices of Eq. (fT4|) for i = J — 1 and 
i = J +1, respectively. On the particle basis, the P matrix for these coupled states is 4x4. 



(15) 



In Ref. |23l], the imaginary parts of the square wells were assumed to be equal for 1 = 
and / = 1 in each partial wave. We take these to be different here, because this choice gives 
a better fit to the more recent high-quality charge-exchange data. The fitted values of the P- 
matrix parameters are given in Table [H The fit to 3749 scattering data requires a total of 46 



TABLE I. P-matrix parameters for the different partial waves. Vq and Vi are the real parts and Wq 
and Wi are the imaginary parts of the short-range sphcrical-well potential, for isospin / = and 
1=1, respectively. The values of the mixing angles 6ij that parametrize the off-diagonal P matrix 
for the partial waves coupled by the tensor force are: ^oi = 7.6° it 0.4° and 6ii = —10.7° it 0.8° for 
the ^Si-^Di waves; 602 = 0.0° and 612 = -8.8° ± 1.6° for the ^P2-^F2 waves; 603 = -7.4° ± 0.4° 
and = —6.9° ± 1.4° for the ^D^-^Gs waves. The quoted errors are defined as the change in each 
parameter that gives a rise in Xmin 1 when the remaining parameters are refitted. 



Partial wave 


Vo (MeV) 


Wo (MeV) 


Vi (MeV) 


Wi (MeV) 







-161.7(25.2) 


-516.1(19.4) 


-132.8(19.9) 




-135.6(9.5) 


-166.9(8.3) 


33.6(5.7) 


-166.3(8.0) 


'Pi 





-374.5(29.6) 





-413.8(40.7) 


'Po 


-114.9(10.1) 


-142.8(9.3) 


-164.1(4.5) 


-71.9(6.9) 


'Pi 


-78.0(4.2) 


-62.2(3.7) 





-382.2(27.6) 


'P2 


-114.6(5.7) 


-201.4(5.1) 


-41.4(3.0) 


-135.6(5.4) 


'D2 


-277.8(16.2) 


-330.8(27.0) 


-319.6(30.4) 


-482.8(45.8) 


'D, 





-96.6(15.5) 





-129.4(19.7) 


'D2 


-120.7(17.6) 


-95.5(16.8) 





-338.6(27.3) 


'Ds 


-235.7(7.7) 


-181.1(8.4) 


-102.0(9.1) 


-66.6(7.8) 




-510.0(22.9) 


-312.4(35.6) 





-335.3(82.0) 







-356.0(56.6) 


-554.0(26.5) 


-317.1(27.0) 




-498.4(61.0) 


-423.2(46.6) 









P-matrix parameters. Almost all the short-range square-well potentials are attractive. The 
quoted errors reflect the sensitivity of the fit to variations in the corresponding parameters. 
These errors are defined as the change in each parameter that gives a rise in Xnim 1 when 
the remaining parameters are refitted. The lower partial waves all require parameters to 
obtain a good fit. To decide which parameters to keep in the fit, a three-sigma criterion is 
used: When the error turns out to be more than one third of the parameter value, it implies 
that xLin ^ises by less than 9 when the remaining parameters are refitted. In that case the 
parameter is set to zero, i.e. it is left out. Because of the centrifugal barrier, the fit becomes 
progressively less sensitive to short-range parameters for the higher-^ partial waves. We 
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assume the parameters in these partial waves to be equal to the ones in similar lower partial 
waves. For example, the parameters for the ^F^ and ^G^ waves are taken to be the same as 
the ones for ^D2; the ones for and ^i^s are the same as the ones for ^^3; and the ones 
for ^G^ and ^Hq are the same as the ones for ^^4. We include the partial waves as high as 
J = 12, which is for instance needed to describe the forward "spike" in the charge-exchange 
differential cross section. 



IV. THE LONG-RANGE ANTINUCLEON-NUCLEON POTENTIAL 



The potential tail for r > b includes the electromagnetic and the strong (nuclear) inter- 
action Vat, where the electromagnetic interaction is the one-photon exchange potential, i.e. 
the Coulomb potential and the magnetic-moment interaction 32| . 



V = Vc + Vmm + Vn . (16) 

In contrast to the NN PWAs, we do not include the vacuum-polarization potential, because 
its effects are negligible, except for very low energies , where there are no pp scattering 



data available. Two-photon exchange effects 36| are not taken into account either. 



The Coulomb potential acts only in the pp channel and is given by the expression 

Voir) = , (17) 

where a' takes care of the main relativistic corrections to the Coulomb potential. It is 
defined by the relativistic Coulomb factor rj = a'Mp/{2p). The magnetic- moment potential 
in the pp channel is given by 

where fip = 1 + Hp = 2.793, with Hp the anomalous magnetic moment of the proton; the 
tensor operator S12 = 3 (Ti ■ r (T2 ■ r — (Ti ■ 0-2, with <Ti and (T2 the spin operators of the 
two nucleons, L is the angular momentum vector, and S = {cri + (T2)/2 the total spin. 
The spin-orbit potential is due to the interaction of the magnetic moment of one particle 
with the charge of the other particle and includes a relativistic correction from the Thomas 
precession. The tensor force is due to the interaction between the magnetic moments of 
the two particles. The magnetic-moment interaction in the nn channel contains only the 
tensor-force part of Eq. (fTSl) with fj^n = i^n = —1.913 and M„. 



The nuclear potential Vn contains the OPE and TPE potentials for AW scattering. Since 
the strong interaction is invariant under charge conjugation C, the AW potential can be 
obtained from the AW potential by using the operator C. If one assumes that isospin 
symmetry SU{2,I) is exact, one can also use the G-parity operator, which is defined as 
G = Cexp(z7r/2), and thus contains charge conjugation and a rotation in isospin space. 
The OPE potential is isospin dependent, while the TPE potential contains both isosp in- 
independent and isospin-dependent parts. When we define the nuclear potential in isospin 
space for the WW system by 

Vn{NN) = W^n-T2 + + W2^n-f2 , (19) 

the potential for the WW system is given by 

Vn{NN) = -W^f^.f2 + V2. + W^2. ri-r2 , (20) 

which implies for elastic and charge-exchange scattering, respectively, 

V^{VV ^nn) = 2 {W^ - W^.) , (21) 

where the factor 2 is due to isospin symmetry. 

The TPE potential for WW scattering has been derived from the effective nonlinear chiral 



Lagrangian density, w. 
symmetry of QCD [28 



lich implements the spontaneously broken SU (2, L) ® SU{2, R) chiral 



37 



. The leading order of this effective Lagrangian density is the 



nonlinear Weinberg model, 

£(0) ^ J^^pM + M + z^7A757^ r • D^] W , (22) 

with the chiral-covariant derivative 

VN = (d'' + ^coT-7rxD'')N , (23) 

where -D^ = D^^d^n/Fj, and D = 1 + t^'^/F^; M is the mass of the nucleon, = 1.269 
is the Gamow- Teller coupling constant in neutron /3 decay, and = 185 MeV is the pion 
decay constant. The subleading-order chiral Lagrangian density is 

£(1) = -N[8ciD-W^^yF^ + 4c3D^-D'' + 2c4a^,r-D''xD'']N . (24) 
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The constant cq = 1 multiplying the Weinberg- Tomozawa NNnn interaction is fixed by chiral 
symmetry. However, the coupling constants Cj {j = 1,3,4) are low-energy constants that 
have to be determined from experimental data. These constants are of order 0{1/M) and 
their values contain contributions from the "integrated-out" heavy hadrons, in particular 
the N- and A-isobars, and the two-pion resonances £(760) and g{770). (The constant C2 
does not contribute to NN scattering at this order.) 



FIG. 1. The Feynman diagrams for one- and two-pion exchange. 

The Feynman diagrams for the OPE and TPE processes are shown in Fig. [H where 
the OPE diagram and the planar- and crossed-box TPE diagrams come from Eq. fl22|) . 
the "triangle" and "football" TPE diagrams containing the Weinberg- Tomozawa NNtttt 
interaction also come from Eq. ( l22i) . while the other "triangle" TPE diagrams containing 
the Cj {j = 1,3,4) NNtth interactions come from Eq. The pion-exchange potentials of 
Eq. (fT9|) contain isospin-independent and isospin-dependent central, spin-spin, tensor, and 
spin-orbit terms. 



Vn = Vc + Vs (Ti- (72 + Vt Si2 + Vso L ■ S 

+ {Wc + Ws(Ti-(T2 + Wt Si2 + Wso L- S) TrT2 



(25) 



where for OPE only the coefficients Ws and Wt are nonzero, and TPE contains in leading 
order only the terms Vs-, Vr, and Wc-, whereas in subleading order all the terms are nonzero. 
The coefficients in Eq. ( I25l) are written in terms of dimensionless functions as 



(26) 



with n = 1 for OPE and n = 2 for TPE, i = C, S, T, SO, and x = m-j^r. We use the 
conventional rationalized "pseudovector" NNtt coupling constant /, normalized such that 



/ 



0.075 



39 



40| . This means that we introduced the scaling mass m^, chosen to be 
numerically equal to the charged-pion mass = 171^^+, and we defined ^ = 111^^/771^. If the 
Goldberger-Treiman relation were exact, one would have that Qa/Ft, = y/Aiif/ms. 
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The OPE potential contains isospin-dependent spin-spin and tensor parts, with 

ws{x) = /?)X , 

wt{x) = {l + x + x^/^)e-''/x^ . (27) 

For the leading- and subleading-order TPE potential, the isospin-independent and the 
isospin-dependent parts can be written as 

Vi{x) = -fi,i(x) + ^fj,2(a;) , 

71 Mp 
2 TTl 

Wi{x) = -Wi^i{x) + -^Wi^2{x) , (28) 

71 Mp 

where the subscript 1 indicates leading-order and the subscript 2 subleading order. The 
leading-order, static TPE potential contains isospin-independent spin-spin and tensor terms 
and an isospin-dependent central term, with 

vs,iix) = l2Ko{2x)/x^ + (12 + 8x^)Ki{2x) /x^ , 
vt,i{x) = -12Kq{2x)/x^ - {Ih + Ax^)Ki{2x) / x'^ , 
wc,i{x) = [cl + lOco - 23 - Ax^) Kq{2x)/x^ 

+ [cl + lOco - 23 + (4£o - I2)x^] Ki{2x)/x^ , (29) 

where cq = CQ/g\ with qa = F-^^f^f/ms and Kn{2x) {n = 0,1) are the modified Bessel 
functions (the hyperbolic Bessel functions) of the second kind, which have asymptotic be- 



havior Kn{2x) ~ ^TxjAxe "^^ for x — ?► oo. The subleading-order potential contains nonstatic 
terms from Eq. (!22|) and the leading-order terms from Eq. (12^ . which can be written as 

^^,2(^) = ELl«fce-27x^ (30) 

and similarly for the Wi^2{x) terms. The coefficients at are listed in Table HIl where we 
defined Cj = CjMp/g\ (j = 1, 3, 4) and Co4 = Cq + 4c4. 

The OPE and TPE potentials for pp — > pp and for pp — nn are now given by Eq. (pTj) . 
In the OPE potential, we take for pp and nn elastic scattering to be the neutral-pion 



mass m^o_and 
of Refs. 



22 



br charge-exchange scattering the charged-pion mass m^+. In the PWAs 



23|, the pion-nucleon coupling constant = fpmT+ fnp-K- was determined 



from the charge-exchange data. In Ref. 22| = 0.0751(17) was found, and in Ref. [23|] 



= 0.0732(11). The values were consistent with the values for f^^^o and found in the 
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TABLE II. The coefficients of the sub 
spin, tensor, and spin-orbit terms 
C04 = Co + 4C4. 



fiading-order TPE potential of Eq. (|30|) for the central, spin- 
28| |: we define cq = co/g\] cj = CjMp/g\ for j = 1,3,4, and 





ai 


0-2 


^3 


04 


05 


06 




3/4 


9 + 48ci + 24c3 


27 + 96ci + 96c3 


99/2 + 48C1 + 240C3 


54 + 288c3 


27 + 144c3 






-3 


-9 


-33/2 


-18 


-9 


VT,2 




3/2 


27/4 


15 


18 


9 


VSO,2 






-12 


-36 


-48 


-24 


WC,2 


3/2 


4-2C0 


14 - 8co 


31 - 20co 


36 - 24co 


18 - 12co 


WS,2 




-2/3 


-14/3 + 8C04/3 


-31/3 + 20C04/3 


-12 + 8C04 


-6 + 4co4 


WT,2 




1/3 


17/6 - 4co4/3 


26/3 - I6C04/3 


12 - 8co4 


6 - 4co4 


WS0,2 








8-8co 


16 - 16co 


8-8co 



•pj) and n-p PWAs [39|, resulting in the recommended value p = /^^^ = 0.0750(9) for the 
pion-nucleon coupling constant, with no significant evidence for isospin breaking 40|. We 
have taken here the values fpp^o = 0.075 and = 0.075 for the OPE potential for elastic and 
charge-exchange scattering, respectively. In the TPE potential we use for the average 
pion mass {2mj^+ + m^o)/3 = 138.04 MeV and the charge-independent coupling constant 
= Inntt — 0.075. The strong potentials for nn — )■ nn and nn — pp are equal to the ones 
for pp — )■ pp and pp — )■ nn, respectively. 

The values of Cj {j = 1,3,4) were determined in the pp and np PWAs 28|, |29|. The 
Ci term in Eq. f l24p breaks chiral symmetry explicitly, since it is proportional to m^. The 



value of Ci cannot be determined accurately from the NN data. It was fixed theoretically 
at Ci = — 0.76/GeV by assuming a value for the pion-nucleon sigma term 28|. We take 
the same value here. It is interesting, however, to probe the sensitivity of our results to 
variations in C3 and C4. It is difficult to determine C3 and C4 and their statistical errors by 
a fit to the database. Since they are parameters in the long-range interaction for r > b, 
this would require that for each small step in varying C3 or C4, the Schrodinger equation 
would have to be solved for all the energies. However, we found that very good results 
were obtained for the values C3 = — 5.8/GeV and C4 = 4.0/GeV, where we estimate the 
uncertainties to be of the order of 0.5. This means that the values we found are remarkably 

13 



consistent with the values determined in the pp PWA to 350 MeV: C3 = — 5.08(28)/GeV and 
C4 = 4.70(70)/GeV 13] . In the pp and np PWA to 500 MeV the values C3 = -4.78(10)/GeV 
and C4 = 3.96(22)/GeV were found 



29|. One could interpret this 
charge conjugation invariance of the TPE interaction. We leave a more careful study of the 
chiral OPE and TPE potential tail in AW scattering for the future. 



The resulting long-range OPE and TPE potentials should be compared to the ones of 



Ref. 



23| where the charge-conjugated version of the high-quality soft-core Nijmegen one- 



boson exchange (OBE) potential [30|, |3l| was used as long-range interaction. In both cases, 
OPE is included, so one should compare TPE to the exchange of the heavy bosons, in 
particular the two-pion resonances £:(760) and g{770). Since the vector mesons have negative 
charge parity, the coupling constants of ^(770) and u;(782) change sign when going from 
nucleons to antinucleons. When we write schematically for the pp potential 

V{pp ^pp) = W^ + Vs + W, + V^ + ... , (31) 

we obtain for the OBE potential for elastic pp — )■ pp and charge-exchange pp — )■ nn scattering 

V{pp -^pp) = W^ + Ve-W,-V^ + ... , 

V{pp ^nn) = 2{W^-W, + ...) , (32) 

respectively. This should be compared to Eq. (^1^. It implies that for the AW case the 
central potential is relatively weak, because there is a cancellation between the repulsion due 
to the vector mesons and the attraction due to the scalar mesons, there is a strong coherent 
spin-orbit force from the exchange of the scalar and vector mesons, and the tensor forces 
due to OPE and g{770) exchange have opposite sign. For the WW case, a strong coherent 
central attraction results due to scalar- and vector-meson exchange and a relatively weak 
spin-orbit potential. Moreover, a strong coherent tensor potential acts in WW due to OPE 
and ^^(770) exchange. This strong tensor force dominates the charge-exchange pp — )■ nn and 



strangeness-exchange pp — )■ AA processes, where no neutral mesons can be exchanged 



41 



4^. 



The chiral TPE potential in subleading order has qualitatively a number of similar fea- 
tures. Because the values of C3 and C4 are large, the corresponding "triangle" diagrams with 
an NNtttt "seagull" interaction lead to relatively strong potentials. The C3 term gives rise 
to a strong central attraction, while the C4 term gives rise to a strong tensor force with the 

14 



same sign as the tensor force due to OPE. This results in a strong attractive central force 
in the elastic process pp — )■ pp and a strong coherent tensor force in the charge-exchange 
process pp fin. This can be understood because the C3 and C4 terms contain effects 
from "integrated-out" £(760) scalar-isoscalar and f?(770) vector-isovector mesons, respec- 
tively. In fact, these two mesons are prominent broad two-pion resonances. In the potential 
of Refs. 



30 



31| their widths are treated in a two-pole approximation; the lowest-mass poles 



correspond to mesons of masses of about 550 and 650 MeV, respectively, resulting in rela- 
tively long-range potentials. 



V. ANTIPROTON-PROTON DATABASE AND STATISTICS 



The antiproton-proton database was constructed for the first time in Ref. 23|. It included 
all available scattering data below antiproton laboratory momentum 925 MeV/c published 
up to early 1993 in a regular physics journal, i.e. total and annihilation cross sections, differ- 
ential cross sections and analyzing powers for elastic and charge-exchange scattering, total 
cross sections for charge-exchange scattering, and (very few) differential depolarizations for 
elastic scattering. At that time, most of the experiments at LEAR were finished. However, 
some more data sets were published after the completion of the PWA of Ref. 23|. We in- 
clude these data sets here, along with a few data sets for which the numerical values were 
not available back then. The present, new database is summarized in Table IIIII The data 



sets that were not included in Ref. 



23| are marked with an asterisk in the left column. We 



always consult the original publications for information about the data and their statistical 
and systematic uncertainties. 

Statistical tools are an essential part of the data analysis in a PWA. We use exactly the 
same methods as in the NN PWAs 25|]. We mention here only the main relevant points, more 



details can be found in Ref. 



23| . We perform a least-squares fit of the model parameters to 



the total database, which contains individual data sets labelled by A. One data set contains 
A^^ individual data points labelled by i. The of the fit is correspondingly defined as 

2" 



X 



i=l 



(33) 



where p is the parameter vector with A'par entries, MA,i{p) is the value predicted by the 
model for the measured observable Eai labelled i in set A with statistical error eAi (in 
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several cases, point-to-point systematic errors were added in quadrature to the statistical 
errors in the experimental papers). In most cases, the data sets have an overall normalization 
uncertainty, denoted by eA,o, specified by the experimentahsts. For each of these sets we 
introduce a normalization parameter ua that multiplies the measured values E^^i of the 
entire set. In the case that the experimental data sets are only relative, or in the case that 
the normalization error was underestimated, the error eA,o is taken to be oo (in practice very 
large) and the corresponding normalization parameter ua is "fioated." The contributions to 

of these normalizations are then zero. In a few cases the normalizations are absolute, i.e. 
eA,o = 0, and the contributions to of these normalizations are again zero. 

By using a sophisticated numerical fitting code, the value of x^ip) is minimized with 
respect to the model parameters. By using the definition Eq. fl55]) . the normalization pa- 
rameters are adjusted implicitly. According to the theory of least-squares fitting, the expec- 
tation value of the minimum is (Xmin) — ^df ± V^^df, where N^^f is the number of degrees 
of freedom, provided the data points are distributed statistically {i.e. they do not contain 
systematic errors) and provided they are Gaussian (which is the case for counting experi- 
ments with enough events per bin). The error matrix E of the model parameters is defined 



where Pmin the values of the model parameters in the minimum value of x^- The error 
matrix allows us to determine the error in the model parameter as y/E^a- This error 
corresponds to the variation in that parameter that gives a rise in Xmm of 1 when the 
remaining parameters are refitted. As mentioned in Sec. Illlt when the error is more than 
one third of the parameter value, it implies that Xmm rises by less than 9 when the remaining 
parameters are refitted. In that case the parameter is set to zero, i.e. it is left out. The 
error matrix allows us also to provide statistical uncertainties on our predictions for the 
observables. In the plots of the differential observables below, the PWA result is given as 
a full red line with an area bordered by blue dotted lines that indicate the one-standard- 
deviation uncertainty in the prediction. 



by 




(34) 
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TABLE III. Reference table of antiproton-proton scattering data with piab < 923 MeV/c. The 



23[, because 



asterisks in the leftmost column indicate the data sets that were not included in Ref. 
the data are more recent or because the values of the data points were not available. The meanings 
of the superscripts in the heading and the comments in the rightmost column are given at the end 
of the table. 



J^lab 


No.'^ 




Norm 


Pred. 








(MeV/c) 


type'' 


Amin 


error^ 


norm'^ 


Rejected*^ 


Ref. 


Comment 


119.0-923.0 


50 fJcc 


46.5 


4% 


1.058 


<385.0, #=8; 468.0 


[45] 


k, m 


176.8-396.1 


5 C^ann 


9.4 


4.4% 


0.949 


176.8 


[46] 




181.0 


46 do-el 




5% 




all 


[47,48] 


j, 1, o 


183.0 


13 dfJce 


13.3 


5% 


1.002 


0.940,-0.170,-0.574 


[49] 




194.8 


19 do-el 




4% 




all 


m 


f, i, o 


200.0-588.2 


48 0-ann 


52.5 


2.2% 


0.989 




[46,5]J 




221.9-413.2 


45 0-tot 


55.3 


oo 


0.961 


221.9, 229.6, 254.9, 260.8, 
280.3, 289.1, 394.2, norm 


[52] 




233.0 


54 do-el 




5% 




all 


[53] 


f,j 


239.2 


20 do-el 


16.0 


4% 


1.077 


-0.950 


[50] 


o 


272.0 


65 do-el 


61.8 


5% 


1.005 




[53] 


j 


276.0-922.0 


21 0-ce 


26.2 


7.5% 


1.098 




m 


m 


276.9 


20 do-el 


20.9 


4% 


1.027 




m 





287.0 


54 do-el 




5% 




all 


[47,48] 


j, 1, o 


287.0 


14 do-ce 


29.6 


5% 


1.144 




[49] 




310.4 


20 do-el 


30.6 


4% 


1.024 




[50] 


o 


340.9 


20 do-el 


23.3 


4% 


1.033 


-0.950, -0.850 


[50] 


o 


348.7 


38 do-el 


40.9 


4% 


0.973 




[55] 


i, o 


353.3 


119 do-el 


117.6 


5% 


1.007 


0.366 


[56] 


j, o 


355.0-923.0 


36 0-tot 




1.5% 




all 


m 


e, m 


359.0-652.0* 


11 do-el 




2% 




all 


[58] 


t 


369.1 


19 do-el 


16.0 


4% 


1.015 


0.550 


m 


i, o 


374.0 


39 do-el 


27.8 


5% 


1.040 




m 


o 


388.0-598.6 


29 fJtot 


35.2 


oo 


0.964 


504.8, norm 


[60] 




392.4 


19 ddel 




5% 




all 


[61] 


1 
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TABLE III. (Continued). 



Til 1 

Plah 


No. a 






J- 1 C-U.. 








(MpV lr\ 


type 


A,min 


Cl 1 Ui 


i±U± ill 




iVCi . 


V_yUiiiiiiCiiL 






8 fi 
o.u 




1 . luO 


Q^^^ 01^4 944 OQQ 

Kj.aOO^ U.c/tJ4:, VJ.Z(4:^, U.Uc/e/, 

O^Q 9*^0 fil 


[fi9l 


W 




98 


1 Q 4 


n 8% 

u.o /o 


1 007 


7*^7 4 








OO tic el 


ou.o 


4% 


Q74 

U.y / 4: 




|DOJ 


1, 






9Q Q 


OO 


7^7 


liUi 111 


ffi5l 


n 




1 1 Q rln- , 


1 no 1^ 


d /O 


1 009 


QQl 7^0 ^7Q 

U.yjyl, U.iOU, U.Oic/ 


|DOJ 




41 1 9 


OO LLUqI 


oo.u 


fj /O 


QQ8 
u.yyo 


_0 87^^ —0 Q9^ 


f59l 

|oyj 


1, 


41 *^ 4 


7 rlrr , 


4 ^ 


o /o 


1 04*^ 


QQ9 


|00J 


J; 




7 rlrr , 




o /o 




nil 

dii 


|DDJ 


6) J) 


49R n 


iU aCTce 


Q fi 

y .0 


90% 


1 1 70 




ffi7l 






7 rlrr i 
' UCTel 


1 9 


/o 


1 007 


QQ9 


|DOJ 


j) 




97 r\rr i 

^ / LLC el 




10% 




all 

dii 


ffi8l 


1 

i 




94 4 1 


ou.u 


OO 


1 'i7Q 
< y 




ffi8l 






Oc/ LLC el 


4n 8 


o /o 


1 OOfi 




f59l 
[oyj 





440 R 


■^8 r]rr i 


48 7 


o /o 


1 094 


79"^ 


f59l 
|oyj 


1, 


AAA 1 


■^8 r\rr i 

OO LLC el 


48 9 


4% 


qfi7 


—0 ^7"^ 


f55l 

|00J 


1, 


44fi n 


1 1 Q rln- , 


1 1 1^ 7 


/o 


QQ8 
u.yyo 




|DDJ 


j) 


447 1 




U. 1 


d /O 


1 0*^8 
1 .uoo 


QQ9 
u.y 


|00J 






8 rlrr , 


9 9 


/o 


Q8fi 
u.you 


QQfi 
u.y c/u 


|DDJ 


j) 


4f^7 ^ 


Oct U.CJel 


Ol.O 


4% 
'I /o 


1 01 Q 


—0 Q9'=l 


f^5l 


1, 


4fi7 8 


Hit 1 


23 5 


5% 


1 033 




f59l 


\J 


469.2 


8 dcJei 


8.0 


5% 


1.004 


0.996 


[66] 


j, 


A7Q 


1 1 Q rln- 1 


1 HQ 
±uy .(J 


fj /o 


Q89 


Q1 Q 87"^ fiQ7 
u.yxy, u.o/o, u.uy/ 


f56l 


J) 


480 n 


1 n r\rr 


1 n 1 


OO 


1 1 1 
± . ± xo 




|oyj 


(T 

8 


/LSI 9 




D.O 


f^% 

/O 


1 0'?7 


QQfi 
U.yyo 


[DOJ 




4Qn 1 


■^7 rlrr i 
o 1 U-U el 




5% 




all 

dii 


[fill 


1 

1 


400 1 


1 '\ r\n- 


lO. 1 


/o(^o.y /o J 


1 001 


QQ9 1 '^Rl 
—0 "^fifi 

U.fJUU 


ffi9l 

|DZJ 


w 


490.6 


39 drr i 


46.3 


5% 

tJ /U 


0.963 




f59l 


Q 


492.7 


8 dcTel 


3.9 


5% 


1.003 


0.996 


[66] 


j, 


497.0 


14 ^y,el 


10.3 


4.5% 


1.004 




[70, UJ 




498.7 


37 dcJei 


28.4 


4% 


0.989 




[55] 


i, 


503.8 


8 dcTel 


12.6 


5% 


1.034 


0.996 


[66] 


j, 


504.7 


39 dcJei 




5% 




all 


[59] 


e, 
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TABLE III. (Continued). 



Plab 


No.^ 




Norm 


Pred. 








(MeV/c) 


type^ 




error'^ 


norm^ 


Rejected"^ 


Ref. 


Comment 


505.0 


54 dcTel 




5% 




all 


[47,48] 


j, 1, o 


505.0 


14 do-ce 


17.3 


5% 


1.021 




[49j 




508.0 


119 dcJei 


106.2 


5% 


0.998 


0.663, 0.530 


[56j 


j, o 


508.9 


39 dfJei 


29.6 


5% 


1.005 




[59j 


o 


516.0 


8 dad 


5.1 


5% 


1.006 


0.996 


[66j 


j, o 


523.0 


15 AyQl 


11.5 


4.5% 


1.037 




[70, 7JJ 




524.8 


36 dfJei 


32.2 


4% 


1.004 




[55J 


i, o 


525.9 


39 do"ei 


42.5 


5% 


1.033 




[59] 


o 


528.2 


8 dcJel 


2.9 


5% 


0.993 


0.996 


[66] 


j, o 


533.6 


119 dcJel 


126.1 


5% 


1.012 


0.892 


[56] 


j, o 


537.0 


10 dfJce 


12.5 


oo 


1.179 




[69] 


g 


540.6 


8 dcJel 


10.7 


5% 


1.004 


0.996 


[66] 


j, o 


543.2 


39 dcTei 


43.9 


5% 


1.051 




[59] 


o 


544.0 


33 dcTel 




10% 




all 


[68] 


1 


544.0 


30 Ay^el 




5% 




all 


[68] 


f, g, o 


546.0* 


12 dcJce 


12.7 


15% 


1.219 




[72] 




546.0* 


2 dcTce 


1.0 


15% 


1.024 




[72] 




546.0 


23 Ay^QQ 


23.3 


4% 


0.966 


-0.250 


[73] 




546.0* 


13 AyQQ 




4% 




all 


[7^ 


f 


546.0* 


7 D 


4.9 


— 


— 




[7^ 




549.4 


10 dcJce 


5.9 


20% 


1.219 




[62 




550.0 


67 ddei 


80.2 


5% 


0.978 


0.997, 0.996, 0.995 
0.910, 0.883 


[75] 


j 


553.1 


34 dcTei 


37.1 


4% 


0.967 




[55] 


i, o 


553.4 


8 dcJci 


2.2 


5% 


1.008 


0.996, 0.972 


[66] 


j, o 


556.9 


119 dcTci 


124.1 


5% 


1.002 


0.908 


[56] 


j, o 


558.5 


39 do"ci 


45.2 


5% 


1.021 




[59] 


o 
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TABLE III. (Continued). 





No.^ 




Norm 


Pred. 








(MeV/c) 


tvpe^ 




error^ 


norm'^ 


Rejected'^ 


Ref. 


Comment 


565.5 


8 dfjoi 


5.4 


5% 


0.994 


0.996 


[661 


i. o 


568.4 


37 dcjgi 


35.4 


5% 


1.025 


-0.675, -0.825 


[591 


i. o 


577.2 


36 dcTgi 


33.9 


4% 


0.969 




[551 


i, o 


578.1 


9 dcToi 


6.0 


5% 


1.001 


0.999 


[661 


i. o 


578.3 


119 dcJci 


133.0 


5% 


1.022 




[561 


i. o 


584.0 


10 dcJce 


6.3 


oo 


1.043 




[691 


O 


590.0 


39 dcTgi 




5% 




all 


[47, 481 


i. 1. o 


590.0 


15 dcTr-o 


23.4 


5% 


1.030 


0.996 


[491 




591.2 


9 dcToi 


6.1 


5% 


1.016 


0.999 


[661 


i. o 


591.2 


39 dcJci 




5% 




all 


[611 


1 


591.2 


15 dcTr-r, 


11.3 


5%(7.8%) 


1.030 


-0.358, -0.545 


[621 


w 


596.5 


38 dcTgi 


46.0 


5% 


1.059 




[591 


o 


599.2 


33 dcTgi 


15.8 


4% 


0.983 




[55.1 


i. o 


601.0* 


47 dcJro 


47.5 


3% 


1.035 




[761 




601.5* 


47 dcJcc 


37.8 


3% 


1.074 




[771 




604.0 


9 dcToi 


8.3 


5% 


0.975 


0.999 


[661 


i. o 


615.0 


38 dcTgi 


55.0 


5% 


1.036 




[591 


o 


617.0 


9 dcToi 


8.0 


5% 


0.944 


0.998 


[661 


i. o 


630.0 


10 dcTr-o 


7.1 


oo 


1.046 




[691 


E 
o 


630.9 


9 dcei 


4.5 


5% 


0.991 


0.999 


[661 


i, o 


639.6 


38 duei 




5% 




all 


[591 


e, o 


644.7 


9 dcToi 


9.2 


5% 


0.961 


0.998 


[661 


i. o 


656.0* 


10 dCTr-o 


12.9 


15% 


1.220 




[721 




656.0* 


7 dcJce 


14.6 


15% 


1.165 




[72] 




656.0 


17 AyQQ 


11.2 


4% 


0.982 




[78] 




656.0 


21 ^y^ce 


23.5 


4% 


0.956 




[73] 




658.1 


38 ddei 


44.6 


5% 


0.963 


0.225, -0.675 


[59] 


o 
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TABLE III. (Continued). 



Plsb 


No.^ 




Norm 


Pred. 








(MeV/c) 


type*" 


■ 


error"^ 


norm^ 


Rejected"^ 


Ref. 


Comment 


658.6 


9 dfjci 


8.9 


5% 


0.999 


0.998 


[66] 


j, o 


670.0 


10 dcJce 


6.2 


oo 


1.150 




[69] 


g 


671.5 


9 dcJei 


3.8 


5% 


0.981 


0.998 


[66] 


j, o 


679.0 


26 dcJei 




oo 


• 


all 


m 


h, 1 


679.0 


27 ^y,el 


25.1 


4.5% 


1.005 


0.540 


[70, 71J 




679.0 


1 Dyyel 


3.2 


— 


— 




[79] 


q 


679.1 




4.1 


5% 


0.984 




[80] 


o 


680.1 


38 dcJci 


39.2 


5% 


0.990 




[59] 


o 


686.1 


9 dcTel 


4.3 


5% 


0.980 


0.998 


[66] 


j, o 


689.0 


39 dcTei 




5% 




all 


m 


1 


689.0 


16 ddce 


17.9 


5%(7.6%) 


0.961 


0.998, 0.981 


m 


w 


690.0 


89 do-el 


94.5 


4% 


0.978 


0.370 


m 




693.0 


24 dcJce 


37.8 


10% 


1.103 




[78] 


r 


693.0 


17 dcTce 


20.4 


10% 


1.041 




[78] 


r 


696.1 


21 dcJci 


18.5 


4% 


1.016 




[82] 




696.1 


16 ddce 


15.3 


4% 


1.031 




[82] 




697.0 


24 dcTei 




10% 




all 


[83] 


1 


697.0 


33 ^yfil 


20.8 


oo 


1.213 


0.629, norm 


[83] 


o 


698.0 


10 ddce 


7.0 


oo 


1.195 




[69] 


g 


700.0 


4 ^J/,el 


1.3 


5% 


0.997 




m 


o 


701.1 


9 dcJel 


3.8 


5% 


0.994 


0.998 


[66] 


j, o 


715.3 


9 dcJei 


10.6 


5% 


0.997 


0.998 


[66] 


j, o 


728.0 


10 ddce 


3.1 


oo 


1.065 




[69] 


g 


730.0* 


40 ddce 




5% 




all 


[85] 


f, s, u 


757.0 


72 daei 


81.4 


5% 


1.023 


0.999, 0.997, 
















0.996, 0.991 


[75] 


j 


760.0* 


24 do-ce 




5% 




all 


[86] 
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TABLE III. (Continued). 





No.^ 




Norm 


Pred. 








(MeV/c) 




Amin 


6rror^ 






Ref. 


Comment 


767.0* 


10 d(T 


9.1 


15% 


1.441 




f72l 




767.0* 


8 rl(T^^ 


9.6 


15% 


1.231 




[721 




767.0 


22 A, rf 


28.0 


4% 


1.106 




[731 




780.5 


39 dfJel 




5% 




all 


[611 


1 


780.5 


1 Fi cltT^^ 


6.7 


5%f7 1%) 


0.963 


982 —0 143 


[621 


w 


783.0 


30 da 1 




OO 




all 


[711 


h 1 
11, 1 


783.0 


30 yl 1 


36.2 


4 5% 


1.068 


—0 300 —0 420 


[70 711 




783.0 


3 D 1 


6.4 








[791 


n 


790.0 


95 dfx 1 


95.3 


4% 


1.020 




[811 




860.0 


95 dfx 1 

■ / 'J \JLyJ qY 


61.0 


4% 


1.019 


0.510 


[811 




875.0* 


10 dn 


8.1 


15% 


1.366 




[721 




875.0* 


10 dn 


20.9 


15% 


1.183 




[721 




875.0* 


12 dcJce 


8.4 


15% 


1.417 




[871 




875.0 


23 A, 


12.1 


4% 


1.050 




[731 




875.0* 


19 A, 


19.2 


4% 


0.972 




[871 




875.0* 


13 A, r-n 


14.0 


15% 


1.089 




[871 




875.0* 


9 D 

yy 


5.1 








[74, 881 


Q 


875.0* 


5 K.. 

^^yy^cc 


5.9 








[891 


n 


886.0 


34 dcTei 




OO 




all 


[I3J 


h, 1 


886.0 


34 ^y,el 


34.1 


4.5% 


1.023 




[70, UJ 




886.0 




1.5 








[79] 


q 


910.0 


19 dcTel 




OO 




all 


[90] 




910.0 


21 ^j,,el 


12.9 


5% 


0.990 




[90] 





a The number includes all published data, except those given as O.OiO.O (see Comment 
i), and those having piab > 923 MeV/c (see Comment m). 



b The subscripts "el" and "ce" denote observables in the elastic pp — pp and charge- 
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exchange pp — )■ nn reactions, respectively, "da" denotes a differential cross section 
dcx/dQ, "Aj^" a polarization-type datum (asymmetry or analyzing power), ^^Dyy" a 
depolarization type datum, and ''Kyy" a spin-transfer type datum. "o"tot" stands 
for total cross section, "cTann" for total annihilation cross section, and "(Tee" for total 
charge-exchange cross section. 

c Normalization, predicted by the analysis, with which the experimental values should 
be multiplied before comparison with the theoretical values. 

d Tabulated is piab in MeV/c, cos^, "norm" or "all." The notation "<385.0, #=8" e.g. 
means that the 8 points with piab <385.0 MeV/c are rejected. The "norm" means that 
the given normalization is rejected and a "floated" normalization is used instead. The 
"all" means that all of the data points in this set are rejected. 



/2. . 



mm 



e Group rejected due to improbable low x 

f Group rejected due to improbable high Xmin- 

g "Floated" normalization. Data are relative only. 

h Normalization "floated" by us, since the norm contributes much more than 9 to Xmi 

i Data points given as O.OiO.O not included. 

j Coulomb-nuclear interference measurement. Data points in the extreme forward an- 
gular region are rejected when they contain multiple-scattering effects. 

k Data points at low momenta rejected. 

the database. For detailed 



1 Problematic differential cross sections. Not included in 
explanation, see Sec. VIIIB and Tables II and III of Ref. 



23|. 



m Part of a group of data with points having piab > 923 MeV/c. 

n Elastic differential cross sections as a function of momentum taken at backward angle 
cos^ = -0.994. 

o Normalization error assumed by us, since no clear number is stated in the reference, 
p Depolarization data. Not included in the fit, in view of the large error bars. 
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q Normalization error taken to be zero, in view of the large error bars of these data. 



r Data points taken at the same angles averaged. 



s Data taken from the website |http : / /hepdata . cedar .ac.uk 



t Data not available. 

u The momentum is the average of 700 MeV/c and 760 MeV/c. 

V Normalization errors as used in the fitting, as deduced from the experimental articles; 
when not explicitly given, a reasonable value was assumed by us. 

w The x%{y%) notation means that x% is the overall normalization error and y% is the 
point-to-point systematic error. 

The total Xmm is only a global measure for the quality of the fit. In the Appendix, we 
discuss in more detail the statistical quality of the final database, by examining the final 
distribution of the data points and how it compares to theoretical expectations. 

Data selection is a necessary ingredient of a PWA. In PWAs of large amounts of scattering 
data, a significant minority of the data sets turns out to be inconsistent with the rest of 
the database and with the PWA solution. In these cases, the data sets usually suffer from 
large systematic errors, which cannot be traced and corrected for. Examples in our case 
are the elastic differential cross sections measured at LEAR, which are inconsistent among 
themselves and with earlier measurements, and which in many cases cannot even be fitted 



properly with Legendre polynomials, as discussed at length in Ref. [23] . Including these 
flawed data sets would seriously bias the PWA solution. To decide whether a data set or 
an individual data point is acceptable, we use the standard statistical criteria outlined in 



Ref. 25| and already applied in Ref. 23|. They are generalized three-sigma criteria: Any 
single data point with • > 9 is rejected, as well as any data set with significantly too 



high or too low x% according to the limits given in Ref. 25|]. Some data or data groups are 



rejected because of other reasons, as mentioned also in the Comments column of Table IIIII 
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I I . I , I , I , I . I . I , 

100 200 300 400 500 600 700 800 

Pi3b(MeV/c) 



FIG. 2. (Color online) Total cross sections and total annihilation cross sections as function of 



antiproton laboratory momentum. The PWA fit has for Briickner et al. 

points (Tann and Xmii 



iTtot; for Clough et al. 
points cTtot- 



52.5 for 48 points 

Cannj for Nakamura et al. 



2 

min 



35.2 for 28 points fitot; for Bugg et al. 



m xLn = 9-4 for 4 
19.4 for 27 points 



55.3 for 38 



VI. DESCRIPTION OF THE DATA 



The final pp database contains Nobs = 3636 scattering observables. The details for each 
of the data sets can be found in In Table IIIII We need A^par = 46 model (P-matrix) param- 
eters for an optimal fit. In the fit we must determine at the same time A''^ normalization 
parameters, so the total number of free parameters is A^fp = A^par + A^n- Of the A^n normal- 
ization parameters A'ne have a finite error, while the rest, A^^f = ~ -^ne, is the number 
of "floated" normalizations. In our case, the total number of normalizations is 131, but we 
fixed the normalizations for the five depolarization Dyy and for the one spin transfer Kyy 
measurements, because these data sets have relatively large error bars. Therefore, Nj^ = 125. 
Of these, A^nf = 12 normalizations are "floated," either because the data sets are relative 
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FIG. 3. (Color online) Total charge-exchange cross sections ace as function of antiproton laboratory 



momentum. The PWA fit has for Alston-Garnjost et al. 
et al. ^] Xmin = ^6.5 for 41 points. 



26.2 for 21 points; for Hamilton 



only, or because the normalization errors were underestimated in the experimental papers. 
Thus, the number of normalizations with errors is A^ne = — ^nf = 113. This implies that 
the total number of free parameters is A^fp = A^par + = 171, the total number of data is 
A'dat = A'obs + A^nc = 3749, and the number of degrees of freedom is A'df = Adat — A^fp = 3578. 
The fit results in a minimum value of Xmm ~ 3750.6. Therefore, the minimum per 
datum is Xmin/^dat = 1.000, and the minimum per degree of freedom is Xmin/^df = 1.048. 
When the model is perfect and the database is a perfect statistical ensemble, one expects 
(Xmin/^df) = 1.000 ± 0.024, hence our result for Xmin/^df is only two standard deviations 
too high. The quality of the fit implies in particular that the charge- conjugated chiral OPE 
and TPE potential provides an excellent long-range AW interaction. 

A detailed discussion of most of the data sets can be found in Ref. Here we will 

show the results for a number of important data sets, and in particular address the high- 



quality data sets that were not available in Ref. [23| . The data sets in the figures have been 
multiplied by the predicted normalization factors given in Table UTTl The rejected outliers are 
not plotted in the figures. In case point-to-point systematic errors were added in quadrature 
to the statistical errors, we plot these total errors. 
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300 400 500 600 700 800 900 1000 

Piab(MeV/c) 

FIG. 4. (Color online) Elastic differential cross sections da/cKl at backward angle, cosO = —0.994, 
as function of antiproton laboratory momentum. The PWA result is given by the drawn red line 
and the dotted blue lines indicate the one-sigma uncertainty region. The fit has for Alston-Garnjost 
et al. iSa] Xmin = 29.9 for 30 points. 

In Fig. [2] the total cross sections atot and the total annihilation cross sections (Jann are 
plotted as function of puh-, the antiproton momentum in the laboratory frame. For the 
annihilation cross sections, we introduced two different normalization parameters for the 
data taken with a thin and with a thick target, c/. Table IIIII In Fig. [3] the total charge- 
exchange cross sections o"cc are plotted as function of piab- Unfortunately, there are no good 
data that map out the rise of the cross section above the pp — Tin threshold at piab — 99 
MeV/c. In Fig. H] the elastic differential cross sections da/dQ at backward angles with 
cos^^ = —0.994 are plotted as function of the momentum in the laboratory frame. These 
data are described well, but the normalization of the data set was "floated." At low energies, 
the theoretical uncertainty of the PWA is significantly smaller than the errors of the data 
points. 

The partial-wave cross sections for both elastic and charge-exchange scattering at piab = 
200, 400, 600, and 800 MeV/c are given in Table HVl It is clear that, in contrast to AW 
scattering, many partial waves contribute to AW scattering already at low energies. The 
reason is that the A7V potentials, in particular the central and tensor components, are very 
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TABLE IV. Partial-wave elastic and charge-exchange cross sections, total cross sections, and total 



annihilation cross sections, in mb, for piab = 200, 400, 600, and 800 MeV/c. 







PP- 


pp 






pp 


— )■ nn 




Piab (MeV/c) 


200 


400 


600 


800 


200 


400 


600 


800 




15.7 


7.9 


4.1 


2.1 


0.7 


0.1 








0.9 


2.5 


4.5 


5.6 


0.8 


0.1 








0.1 


0.4 


1.4 


3.1 


0.1 


0.3 


0.1 




'F3 




0.1 


0.2 


0.5 




0.1 


0.1 


0.1 










0.1 






0.1 


0.1 


^Po 


4.9 


5.4 


5.0 


3.5 


1.5 


0.8 


0.1 




^Pl 


1.8 


4.9 


4.0 


3.5 


4.9 


2.9 


0.2 


0.1 




0.1 


0.3 


1.0 


1.5 


0.3 


2.4 


2.5 


1.0 






0.1 


0.1 


0.2 




0.4 


1.1 


1.4 








0.1 


0.1 




0.1 


0.3 


0.5 




66.1 


26.0 


13.2 


8.8 


3.0 


1.0 


0.5 


0.2 




0.3 


0.4 


0.2 


0.1 


0.8 


1.5 


1.1 


0.6 




0.3 


0.4 


0.2 


0.1 


2.0 


2.0 


1.2 


0.7 




0.1 


0.5 


0.8 


1.0 


0.1 


0.5 


0.6 


0.4 


^P2 


7.0 


17.0 


13.9 


9.6 


0.9 


1.4 


0.4 


0.1 




0.1 


0.1 


0.1 




0.1 


0.5 


0.5 


0.5 


^F2 ^3P2 


0.1 


0.1 


0.1 




0.3 


0.8 


0.6 


0.5 








0.1 


0.4 






0.1 


0.1 






1.6 


5.9 


7.0 




0.5 


1.3 


0.6 






0.1 


0.1 






0.2 


0.3 


0.3 






0.1 


0.1 






0.3 


0.5 


0.4 


















0.1 








0.3 


0.8 






0.1 


0.3 


'^F.i -^-^Hi 








0.1 






0.2 


0.2 


3 Tj i.3 IP 








U.i 




n 1 
U.i 


U.Z 


U.o 


3 TT 


















Rest 






0.1 


0.3 




0.1 


0.4 


0.8 


Singlet 


16.7 


10.9 


10.2 


11.3 


1.6 


0.6 


0.4 


0.3 


Triplet 


80.8 


56.9 


45.1 


37.1 


14.0 


15.6 


12.1 


9.2 


Total 


97.5 


67.9 


55.3 


48.4 


15.6 


16.2 


12.5 


9.4 






pp - 


■> all 






pp 


mesons 




piab (MeV/c) 


200 


400 


600 


800 


200 


400 


600 


800 




311.2 


192.6 


149.8 


126.4 


198.1 


108.5 


81.9 


68.6 
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FIG. 5. (Color online) Differential cross sections and analyzing powers for elastic scattering as 
function of angle in the center-of-mass system. The PWA result is given by the drawn red line 
the dotted blue lines indicate the one-sigma uncertaint y reg ion. The fit has for Sakamoto et 



ana 
al. 



2 

min 



39.2 for 38 points da/dVt; for Kunne et al. 



Eisenhandler et al. 
points Ay. 



70 



m X 



2 

min 



Bll Xmin = 94.5 for 88 points da/dQ; for Bertini et al. 



25.1 for 26 points Ay] for 
20.8 for 32 



^ Xmin 



strong. The dominance of the tensor force is seen in particular in the charge-exchange 
pp — )■ nn reaction. For low energies of the final-state nn system the strong tensor force leads 
to large cross sections for the transitions i{nn) = i{pp) — 2, in particular ^Di — )■ ^Si and 
— 7- ^P2- This is similar to the strangeness-exchange reaction pp — )■ AA, where these off- 
diagonal tensor- force transitions due to K[494) and K*{892) exchange dominate the cross 
section in the AA threshold region For these transitions, there is a large overlap 

between the wave functions of the initial pp state and the final nn or AA state [4J| at low 
energy. The contributions from the spin-triplet states are much larger than the contributions 
from the spin-singlet states, especially for pp — )■ nn. The total annihilation cross section is 
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FIG. 6. (Color online) Differential cross sections and analyzing powers for elastic scattering as 
function of angle in the center-of-mass system. The PWA result is given by the drawn red line and 
the dotted blue lines indicate the one-sigma uncertainty region. The fit has for Eisenhandler et 



al. 



8l\ Xmin ~ 95-3 for 95 points da/dQ; for Kunne et al. 



2d, 



7l|] Xmin = 36.2 for 28 points Ay. 
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■a 
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FIG. 7. (Color online) Differential cross sections and analyzing powers for elastic scattering as 
function of angle in the center-of-mass system. The PWA result is given by the drawn red line and 
the dotted blue lines indicate the one-sigma uncertainty region. The fit has for Eisenhandler et 



al. 



81| Xmin ~ 61.0 for 94 points da/dCl; for Kunne et al. 



Id, 



7l|] Xmin = 34.1 for 34 points Ay. 



large, and decreases from a fraction of about 2/3 of the total cross section at piab = 200 
MeV/c to about 1/2 of the total cross section at piab = 800 MeV/c. 

In Figs. O El and [7] the differential cross sections da/dQ and the analyzing powers Ay are 
shown for elastic scattering pp pp at momenta near 690, 790, and 860 MeV/c, respectively. 
In general, the uncertainty on the PWA prediction for the differential cross sections is 
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FIG. 8. (Color online) Differential depolarizations Dyy for elastic scattering as function of angle 
in the center-of-mass system. The PWA result is given by the drawn red line and the dotted blue 
lines indicate the one-sigma uncertainty region. The fit has for Kunne et al. [79] at piab = 679.0 
MeV/c Xmin = 3-2 for 1 point, at piab = 783.0 MeV/c Xmin = 6-4 for 3 points, at puh = 886.0 
MeV/c Xmin ~ ^-^ 1 point. 
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FIG. 9. (Color online) Differential cross sections and analyzing powers for charge-exchange scat- 
tering as function of angle in the center-of-mass system. The PWA result is given by the drawn red 
line and the dotted blue lines indicate the one-sigma uncertainty region. The fit has for Ahmidouch 
et al. [72] Xmin ~ 12-''' 12 points da/dQ at backward angles, Xmin ~ ^-^ ^ points da/dO. at 



forward angles; for Birsa et al. 



73[ y^l^.^ = 23.3 for 22 points Ay. 



determined by the accuracy of the data. For the analyzing powers, on the other hand, 
the theoretical uncertainties are in general smaller than the errors of the data points. The 
theoretical uncertainty is very small for forward angles. For backward angles, where there 
are no data available, this uncertainty increases. Fig. [8] shows the very limited data available 
for the depolarization Dyy for elastic scattering at 679, 783, and 886 MeV/c. There are only 
a few data points in the backward hemisphere and the data points have large error bars. 
In this case, the theoretical uncertainty for the PWA prediction is much smaller than these 
error bars, which implies that there is little new information in these data and that the fit 
would not change significantly if they were left out of the fit. The theoretical uncertainty is 
again very small for forward angles. 

Figs, m [TOl [m and [12] show the differential cross sections da/dVt and the analyzing 
powers Ay for charge-exchange scattering pp — Tin at 546, 656, 767, and 875 MeV/c, 
respectively. Like for the elastic case, one observes that, in general, the uncertainty on 
the PWA prediction for the differential cross sections is determined by the accuracy of the 
data. For the analyzing powers, on the other hand, the theoretical uncertainties are in 
general smaller than the errors of the data points. For some of the differential cross-section 
measurements, we introduced different normalization parameters for the data in the forward 
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COS0 COS0 

FIG. 10. (Color online) Differential cross sections and analyzing powers for charge-exchange scat- 
tering as function of angle in the center-of-mass system. The PWA result is given by the drawn red 
line and the dotted blue lines indicate the one-sigma uncertainty region. The fit has for Ahmidouch 
et al. [72] Xmin ~ 1^-9 10 points da/dQ at backward angles, Xmin ~ ^^-^ points da/dO. at 
forward angles; for Birsa et al. \y\ Xmin ~ ^^■'^ ^^^^ points Ay] for Birsa et al. [73] Xmin ~ ^^-^ 
for 21 points Ay. 



and in the backward hemisphere, which were taken with different detectors. The charge- 
exchange differential cross section is highly anisotropic, because of the contributions of many, 
high-£ partial waves. It has a "spike" at the most forward angles and it is flat at backward 
angles. It exhibits a very typical dip-bump structure at forward angles, which is due to the 



interference of the OPE interaction with a background due to short-range interactions [91]. 
The precise form of this structure evolves rapidly as function of energy, from a rather flat 
plateau structure at 546 MeV/c to a pronounced dip-bump structure at 875 MeV/c. The 
structure was measured accurately at 601 MeV/c by the PS206 experiment at the end of 
the LEAR era 0]. The high-quality charge-exchange differential cross sections from 

[T3l At the time of Ref. {23], only the data at 693 MeV/c shown 



Ref. [77| are shown in Fig. 
in Fig. [13] were available 



78| . but these differential cross sections did not pin down the 



dip-bump structure. The PWA of Ref. [23| predicted a more pronounced structure for this 
data set. 

In Fig. [TUthe few data sets available for the depolarization Dyy at 546 and 875 MeV/c 
and the spin transfer Kyy at 875 MeV/c in charge-exchange scattering are shown. The data 
points have large error bars, and also in this case the theoretical uncertainty for the PWA 
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FIG. 11. (Color online) Differential cross sections and analyzing powers for charge-exchange scat- 
tering as function of angle in the center-of-mass system. The PWA result is given by the drawn red 
line and the dotted blue lines indicate the one-sigma uncertainty region. The fit has for Ahmidouch 
et al. [72] Xmin ~ 9-1 for 10 points da/dQ at backward angles, Xmin ~ 9. 6 for 8 points da/dO. at 



forward angles; for Birsa et al. 



73[ xl,^^ = 28.0 for 22 points Ay. 
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FIG. 12. (Color online) Differential cross sections and analyzing powers for charge-exchange scat- 
tering as function of angle in the center-of-mass system. The PWA result is given by the drawn red 
line and the dotted blue lines indicate the one-sigma uncertainty re gion . The fit has for Lamanna 



et al. 



m x[ 
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min 



— 1 Xmin 



^.4 for 12 points da/dQ; for Ahmidouch et al. 
da/dQ at backward angles, Xmin ~ ^^.9 for 10 points da/dQ at forward angles; for Birsa et al. [7^ 



8.1 for 10 points 

n 



12.1 for 23 points Ay] for Lamanna et al. 



m X 
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min 



19.2 for 19 points Ay at forward 



angles, xl 



14.0 for 13 points Ay at backward angles. 
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Plat, = 693.0 MeV/c 

Birsaetal. (1990) ■ 
Birsaetal. (1990) ■ 
PWA 




FIG. 13. (Color online) Differential cross sections da/dQ for charge-exchange scattering as function 
of angle in the center-of-mass system. The PWA result is given by the drawn red line and the dotted 



n 



blue lines indicate the one-sigma uncertainty region. The fit has for Bressan et al. [77] Xmin ~ ^'^•^ 
for 47 points; for Birsa et al. [78] Xmin ~ ^'''•^ ^4 points at backward angles, Xmin ~ ^0.4 for 17 
points at forward angles. 



prediction is much smaller than these error bars. This demonstrates that spin observables 
are, of course, important, but they improve a good energy-dependent PWA only if they are 
precise enough 2J|. The theoretical uncertainty is again very small for forward angles. 
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Ahmidouch et al. (1995) ' — •- 
PWA 
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-0.5 





cos9 



0.5 



FIG. 14. (Color online) Differential depolarizations and spin transfers for charge-exchange scat- 
tering as function of angle in the center-of-mass system. The PWA result is given by the drawn 
red line and the dotted blue lines indicate the one-sigma uncertainty region. The fit has for Ah 
74l ] at piab = 546.0 MeV/c Xmin ~ ^-^ for 7 points Dyy; for Ahmidouch et al. 



midouch et al. 
(Birsa et al. 




at piab = 875.0 MeV/c x 



2 

min 



5.1 for 9 points Dyy; for Ahmidouch et al. 



74| 



891] at 



Piab = 875.0 MeV/c xl 



5.9 for 5 points Kyy. 



36 



VII. PHASE-SHIFT AND INELASTICITY PARAMETERS 



In this section we present results for the AW S matrix. The S matrix for the coupled pp 
and nn channels from our PWA suffices to construct the complete scattering amplitudes and 
hence the observables. For the uncoupled partial waves with £ = J, s = 0, 1 or £ = 1, J = 0, 
the 5* matrix is 2 x 2, while for the coupled partial waves with i = J± 1 (J > 1), s = 1, it is a 
4x4 matrix. We give numerical values at a number of momenta. Other results are available 
upon request. The matrix elements of S^+tv fo^' different partial waves for the elastic and 
charge-exchange reactions are given in Tables IV l IVll and IVIII for piab = 100 to 1000 MeV/c. 
The S matrices are symmetric for the coupled partial waves in the case of elastic pp and nn 
scattering, but they are not symmetric in the case of charge-exchange scattering pp ^ nn, 
as one can see from Table |Vl] for pp pp and from Table IVIII for pp — )■ nn. 

For illustrative purposes we also present phase-shift and inelasticity parameters assuming 
that isospin symmetry is exact (we take then the average nucleon and pion mass and set the 
electromagnetic interaction to zero). In that case, the parametrization of the S matrix can 
be done in a transparent way, similar to the procedures used for AW scattering (above the 
pion-production threshold). 

For the uncoupled partial waves with ^ = J, s = 0, 1 or £ = 1, J = 0, the S matrix, for 
isospin / = Oor/ = l, isalxl matrix that can be written as 

= r]exp{2i6) , (35) 

where 6 is the phase shift and t] {0 < rj < 1) is the inelasticity due to the annihilation 
into mesonic channels. The 5* matrix for the uncoupled waves is thus given in terms of two 
parameters, which are functions of energy. For high values of £, where there is almost no 
annihilation, r/ — > 1. 

For the partial waves with £=J±1(J>1), s = l, coupled by a tensor force, the 
S*- matrix, for isospin / = or / = 1, is a 2 x 2 matrix that can be parametrized by the 



generalized "bar-phase" convention 92 1 



S"^ = exp{i5) exp{iejax) exp{iejax) exp{i5) , (36) 

where 5 is a 2 x 2 diagonal matrix with real entries and Sj^i^j, and ej is the mixing 

angle for the coupled partial waves; is the first Pauli matrix. The matrix H"^ is used to 
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parametrize the inelasticity. Different ways to write can be found in the hterature. We 
will follow the parametrization of Ref . [93j , in which one writes 

/ r/j_ij \ 

H'' = exp{-iujay) ' exp{iujjay) , (37) 

where ?7j-i,j and are real numbers with < ?7j=pi,j < 1, and Uj is the mixing angle 

for the inelasticity; ay is the second Pauli matrix. The S matrix for these coupled waves is 
thus given in terms of six parameters. 

From the numerical results of the PWA, one has to extract for each energy the phase-shift 
and inelasticity parameters from the numerical values of the S matrix. For the uncoupled 
partial waves this is easy. In order to obtain the phase-shift and inelasticity parameters for 



the coupled partial waves, the algorithm of Ref. [92| is used. One can write the S matrix as 



_ I ^11 exp(2i5ii) ii?i2 exp(2i5i2) ^^^^ 
ii?i2exp(2z5i2) R22 exp{2i622) 

where Rij and Sij are real numbers. When one defines the auxiliary phases 

&a = — , 

Ob = 622 - , (39) 

6' = 611 + 622 — 2(5i2 , 



it follows that 



sin 2(5' 

tan2^„ + 0, - ^2 



R11R22 + RI2 cos 25' 



1-12 ' 

tan(^, - Ob) = 7 tan(^„ + 6^) . (40) 
Ku + -K22 

From this the phase-shift parameters 5j-i,j and can be obtained. The mixing angle 

£j is given by 

- _ 2R,2Cos{e, + d,-S') 

/?iicos2^„ + /?22Cos2^b ■ ^ ^ 

The elements of the matrix H'^ can then be related to the parameters obtained. One finds 

2Hu cos2ej = -Rii(l + cos2ej) cos26'a + -R22(l — cos2e j) cos29b , 
2i722Cos2£j = -cos2£j)cos2^„ + i?22(l + cos2£j)cos2^6 , (42) 

if 12 cos 2ej = R12 sm{6' -9a- Ob) , 
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from which one can determine the values of Hu, H22, and Hi2. By using Eq. ( 137|) . the 
remaining parameters rjj^i^j, rjj^i^j, and uj can be obtained via 

?7j_i,j + ?7j+i,j = TtH-^ , 

rij^ijrij+ij = detH\ (43) 
tan 2 a; J = 2/712/(^11-^22) ■ 

If one extracts the values of the parameters for one single energy, there can be ambigu- 
ities In order to ensure continuity as function of energy one can always change the 
values of these parameters in such a way that the corresponding S'-matrix elements are not 
changed. In the case of uncoupled partial waves, one can change 6 by 180° and keep rj 
unchanged, as can be seen from Eq. 1^5^. In the case of the coupled partial waves, for 
instance, one can change or by 180° and at the same time change the signs of 
ej and coj, while keeping rjj^i^j and rjj+i^j unchanged; one can also change both and 
by 180° at the same time and keep 77j-i,j, "/^j+i.j, £j, and coj unchanged. In the limit 
where rjj^ij = 1, 5j=pi,j = 0, and ej = 0, one can choose cuj = in order to keep continuity, 
although uj can take any value in this case, but the corresponding S-matrix elements are 
unchanged. 

The results of the phase-shift and inelasticity parameters are given in Tables IVIIIl IIXI 
and |X] for piab = 100 to 1000 MeV/c. A convenient way to plot the S matrix, or equivalently 
T = (5 — 1)/^, as function of energy is to use Argand diagrams. In Fig. [T5l Argand diagrams 
are shown for the uncoupled partial waves and in Fig. [16] for the coupled ones assuming 
isospin symmetry. 
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TABLE V. iS-matrix elements of the uncoupled partial waves for pp — ^ pp and pp — > nn. 



Piab(MeV/c) 




100 






200 






300 






400 






500 






1 c 
•^o 





oyo — 


U 


193i 





aOi — 


A 
U 


296i 





1 ^^o 
lOZ — 


A 
U 


300i 





OQO 

Uoz — 


A 
U 


242i 


-0 


Uoo — 


A 
U 


157i 




3 p 





000 — 


U 


023i 





00/ — 


A 
U 


199i 





40U — 


A 
U 


287i 







A 
U 


248i 


— 


nn7 
UU/ — 


A 
U 


145i 




1 P 





Jo i -\- 


u 


023i 





yio -h 


n 
U 


054? 





oUO -h 


u 


015?; 





Doo — 


A 
U 


076i 





OOo — 


n 
U 


181i 


t 


3 p 





yoy — 


U 


030i 





oyo — 


A 

U 


093i 







A 

u 


173i 





Dio — 


A 
U 


238* 





OZ^ — 


A 

U 


253i 


1 n 


1 


UUU -h 


U 


002i 





yyo H- 


A 

U 


020i 





n77 1 

y/ / H- 


A 

u 


050i 





yzo H- 


A 
U 


074i 







A 

u 


072i 


Sin 


3 n 


1 


UUU — 


n 


003i 





yyo — 


A 
U 


025i 





no 1 

yoi — 


A 
U 


050i 





nQQ 

yoo — 


A 
U 


064i 





ooy — 


A 
U 


071i 




1 


UUU + 


U 


OOOi 


1 


nnn 1 
UUU + 


A 
U 


004i 





nnn 1 

yyy + 


A 
U 


015i 





yy/ + 


A 
U 


031i 





non 1 

yoy + 


A 
U 


051i 




3 F 


1 


UUU -|- 


U 


OOOi 


1 


nnn 
uuu — 


A 
U 


006i 





nnn 
yyy — 


A 
U 


021i 





nnc^ 

yyo — 


A 
U 


039i 





no ^ 

yoo — 


A 
U 


052i 




G4 


1 


UUU + 


U 


OOOi 


1 


nnn 1 

uuu + 


A 
U 


OOli 


1 


nnn 1 

uuu + 


A 
U 


005i 


1 


nnn 1 
UUU + 


A 

u 


012? 





nnn 1 

yyy + 


A 
U 


020i 






1 


nnn 1 
UUU + 


U 


OOOi 


1 


nnn 

uuu — 


A 
U 


OOli 


1 


nnn 

uuu — 


A 
U 


008i 





nnn 

yyy — 


A 
U 


018* 





nno 

yyo — 


A 
U 


030i 






-0 


noi 1 
Uzi H- 


U 


090i 





UOo H- 


A 
U 


131i 





1 nn 1 
lUU + 


A 
U 


035i 





n7n 
U /U — 


A 
U 


041i 





nnn 

uuy — 


A 
U 


077i 




3 p 
-H) 


-0 


nnT 1 


U 


002i 


-0 


oon 
ZzU — 


A 
U 


020i 


-0 


Qni 
oUi — 


A 
U 


082i 


-0 


QnQ 
oUo — 


A 
U 


120i 


-0 


z4/ — 


A 
U 


115i 




1 P 





UUU -r 


U 


003i 


— 


UU 1 -r 


U 


093i 





00/1 _L 
UU4 -|- 


U 


099i 





097 _L 
Uz / -|- 


U 


066i 





OAI _L 


n 
U 


029i 


e 
le 

T 


^Pl 





000- 





006i 





026- 





228i 





045- 





347i 





043- 





345i 





038- 





258i 




U 


000+ 





UUU^ 


A 

u 


000 + 







A 

— u 


004 + 





[J06l 


A 

— u 


008 + 





Uo4z 


A 

— u 


010 + 





volt 









000+ 





AAA ■ 


A 


000 - 





046? 


A 

— 


002- 







A 


004- 







A 


023 - 





S27i 




u 


000+ 





AAA ■ 


A 




000 + 





005i 


A 




000 + 





AOO ■ 


A 

— 


002 + 





U4zz 


A 

— 


007 + 





A C ' 

058^ 




'F, 


u 


000+ 





UUU2 


A 

u 


000- 





UUOl 


A 

— u 


001- 





U4UZ 


A 

— u 


003- 





Uoy^ 


A 

— u 


006- 











u 


000+ 





AAA™* 


A 




000 + 





OOli 


A 




000 + 





A AO™ 

008i 


A 




000 + 





A1 A™ 


A 

— 


001 + 





AO 1 

031^ 




'G4 


u 


000+ 





AAA™* 





000- 





AA1 ™* 

OOli 


A 


000- 





A1 O™ 


— 


001- 





034i 


— 


002- 





062^ 


Piab(McV/c) 




600 






700 






800 






900 






1000 






'So 


-0 


042 - 





076i 


-0 


013- 





019i 





036 + 





007i 





085 + 





002? 





123- 





026? 




'Po 


-0 


162- 





024i 


-0 


236 + 





098i 


-0 


245 + 





207i 


-0 


213 + 





294i 


-0 


159 + 





356i 




'Pi 





429- 





271i 





289- 





330i 





151- 





349i 





027- 





331i 


— 


075- 





280i 


T 


'Pi 





443- 





222i 





360- 





167i 





281- 





llOi 





215- 





056i 





170- 





007i 


'D2 





721 + 





032i 





601- 





035i 





482- 





llOi 





366- 





176i 





253- 





223j 


Rh 


'D, 





789- 





081i 





714- 





090i 





648- 





090i 





585 - 





077? 





521 - 





056z 


'F^ 





971 + 





077i 





934 + 





105i 





874 + 





126i 





794 + 





130i 





703 + 





112« 




'F, 





969- 





057i 





944- 





055i 





910- 





048i 





870- 





041i 





826- 





036i 




'G, 





998 + 





030i 





996 + 





042i 





992 + 





058i 





984 + 





077i 





969 + 





lOli 




'G, 





994- 





041i 





990- 





048i 





983- 





052i 





974- 





051i 





962- 





047i 




'So 


-0 


057- 





071i 


-0 


111- 





030i 


-0 


138 + 





035i 


-0 


136 + 





109i 


-0 


104 + 





180i 




'Po 


-0 


171- 





084i 


-0 


100- 





061i 


-0 


043- 





058i 


-0 


006 - 





076i 





008 - 





106i 




'Pi 





039- 





OOli 





027- 





019i 





012- 





026i 


-0 


002 - 





024z 


-0 


012- 





016« 


e 
le 

T 


'Pi 





039- 





125i 





044 + 





014i 





063 + 





127« 





100 + 





202i 





149 + 





242z 


'D, 


-0 


009 + 





078i 


-0 


007 + 





065i 


-0 


007 + 





054i 


-0 


009 + 





045i 


-0 


011 + 





040i 


i& 


'D2 





050- 





362i 





073- 





348i 





084- 





294i 





082- 





213i 





071- 





123i 


'Fs 


-0 


017 + 





065i 


-0 


035 + 





062i 


-0 


063 + 





046i 


-0 


097 + 





016i 


-0 


129- 





024i 




'Fs 


-0 


006- 





204i 


-0 


001- 





258i 





012- 





303i 





032- 





330i 





056- 





337i 




'G, 


-0 


001 + 





04H 


-0 


003 + 





048i 


-0 


005 + 





052i 


-0 


008 + 





053i 


-0 


015 + 





053i 




'G, 


-0 


004- 





095i 


-0 


006- 





129i 


-0 


007- 





163i 


-0 


007- 





197z 


-0 


004- 





230z 
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TABLE VI. <S-matrix elements of the coupled partial waves for pp — )■ pp. 



piab(MeV/c) 




100 






200 






300 






400 






500 






Sn 





514-0 


307i 





207-0 


290i 





039-0 


195i 


-0 


039- 





096i 


-0 


071-0 


015i 


CO 


S22 





998 + 


003i 





974 + 


015i 





926 + 


025i 





865 + 





034i 





798 + 


OUi 


CO 


S12 





004-0 


015i 


-0 


010-0 


058i 


-0 


031-0 


093i 


-0 


056- 





108i 


-0 


077-0 


103i 




Su 





971 + 


02H 





807 + 


089i 





558 + 


llli 





355 + 





lOli 





220 + 


095j 




S22 


1 


000 + 


OOOi 





998 + 


003i 





991 + 


008i 





980 + 





012i 





964 + 


018z 


CO 


S12 





000-0 


003i 





007-0 


022i 





016-0 


038i 





014- 





047i 





004-0 


055i 


CO 


Sn 


1 


000 + 


OOOi 





998 + 


OlOi 





981 + 


055i 





919 + 





146i 





791 + 


247i 


CO 

1 

CO 


S22 


1 


000 + 


OOOi 


1 


000 + 


OOOi 





999 + 


003i 





996 + 





007i 





991 + 


OlOi 


CO 


S12 





000 + 


OOOi 





000-0 


005i 





002-0 


020i 





008- 





037i 





019-0 


046i 




Sn 


1 


000 + 


OOOi 


1 


000 + 


OOli 


1 


000 + 


005i 





998 + 





018i 





990 + 


045i 


CO 


S22 


1 


000 + 


OOOi 


1 


000 + 


000* 


1 


000 + 


OOH 





999 + 





003« 





998 + 


005i 


CO 







()()() + 


000/ 





000 - 


001/ 





000 - 


007/ 





000 - 





01.",/ 





002 - 


026i 


Piab(MeV/c) 




600 






700 






800 






900 






1000 






Sn 


-0 


087 + 


047i 


-0 


107 + 


095i 


-0 


141 + 


141i 


-0 


191 + 





193* 


-0 


244 + 


254i 


CO 


S22 





734 + 


051i 





671 + 


047i 





605 + 


026i 





526- 





016i 





432-0 


077j 


CO 


S12 


-0 


092-0 


078i 


-0 


097-0 


040i 


-0 


095-0 


OOli 


-0 


089 + 





027« 


-0 


082 + 


037i 




Sn 





141 + 


097i 





097 + 


103i 





073 + 


UOi 





059 + 





imi 





051 + 


121i 


CO 


S22 





938 + 


031i 





892 + 


048i 





822 + 


066i 





725 + 





07Ai 





610 + 


063i 


CO 


S12 


-0 


009-0 


060i 


-0 


024-0 


058i 


-0 


042-0 


050i 


-0 


062- 





038i 


-0 


085-0 


027i 


CO 


Sn 





642 + 


315i 





518 + 


354i 





437 + 


369i 





401 + 





357i 





403 + 


315i 


CO 

1 

CO 


S22 





986 + 


014i 





981 + 


021i 





971 + 


033i 





952 + 





055i 





914 + 


087i 


CO 


S12 





022 -0 


045i 





014-0 


041i 


-0 


003-0 


037i 


-0 


026- 





032i 


-0 


055-0 


023z 




Sn 





972 + 


085i 





940 + 


132i 





897 + 


177i 





855 + 





211i 





822 + 


229i 


CO 


S22 





996 + 


008i 





994 + 


012i 





990 + 


015i 





987 + 





019i 





981 + 


026j 


CO 


S12 





005-0 


037i 





009-0 


046i 





012-0 


052i 





013- 





056i 





012-0 


059z 
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TABLE VII. S'-matrix elements of the coupled partial waves for pp — >■ nn. 



Piab(MeV/c) 




1 nn 






200 






OUU 






zinn 

4:UU 






OUU 




CO 


^11 


A 

— U 


noi n 


071 
U / IZ 


A 

— U 


154- 





noi 
uy iz 


A 

— U 


1 Ql 1 A 


ni n-j' 

UlUZ 


A 

— U 


1 79 -1- n 


lUZZ 


A 

— U 


1 9"^ -L- n 
iZO -f- u 


iDOZ 


'S'22 


A 
U 


uuu + u 


UUUZ 


A 

— U 


014 + 





UZOZ 


A 

— U 


n/i _L n 
U4y + u 


uoyz 


A 

— U 


noi _L n 
uyi + u 


1 1 /It 

1 i4Z 


A 

— U 


1 90 ^ n 
izy + u 


i04Z 


CO 


'S'12 


A 
U 


UUo — u 


Uziz 


A 
U 


003 - 





1 A Qi' 


A 

— U 


n/1 n 
U4z — u 


ZZoZ 


A 

— U 


noQ n 

uyo — u 


zoyz 


A 

— U 


1 A« n 
140 — u 


ZoUZ 




S21 


A 
U 


nnn -l n 
uuu -t u 


UUUZ 


A 

— U 


012- 





UyoZ 


A 

— U 


uoo — u 


1 / yz 


A 

— U 


1 ns n 
iUo — u 


zzyz 


A 

— U 


104 — u 


Z4oZ 




Sn 


A 

— U 


nni -1- n 
uui -\- u 


UUUZ 


A 

— U 


077 + 





ni '^'i 

UloZ 


A 

— U 


1 K9 n 
IDZ — u 


UZDZ 


A 

— U 


181 n 

ioi — u 


UOoZ 


A 

— U 


1 7n n 
1 ( u — u 


UOOZ 


CO 


S22 


A 
U 


nnn -u n 
uuu -t u 


UUUZ 


A 

— U 


001 + 





UUoZ 


A 

— U 


nn^; -u n 
UUO -|- u 


ni 

UlOZ 


A 

— U 


ni 1 J- n 
Uii -\- u 


UoZZ 


A 

— U 


ni R 1 n 
uio -f- u 


UOUZ 


(In 

CO 


S12 


A 
U 


nnn _l n 
uuu + u 


UUUZ 


A 
U 


008 - 





U40Z 


A 
U 


no7 n 
U/ ( — u 


uyyz 


A 
U 


nQc^ n 
UoO — u 


iooz 


A 
U 


n9A n 
uzo — u 


iOZZ 




S21 


A 
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TABLE VIII. Phase-shift and inelasticity parameters of the uncoupled partial waves assuming 

isospin symmetry. 6 is given in degrees. 



Vc) 


100 


200 


300 


400 


500 


600 


700 


800 


900 


1000 


rj 


0.61 


0.48 


0.39 


0.31 


0.24 


0.18 


0.13 


0.11 


0.13 


0.16 


S 


179.28 


169.97 


157.96 


145.40 


132.28 


117.71 


99.83 


77.61 


56.50 


41.44 


V 


0.73 


0.51 


0.34 


0.20 


0.09 


0.02 


0.10 


0.18 


0.25 


0.31 


5 


-14.44 


-26.94 


-38.41 


-49.03 


-57.58 


-4.58 


4.39 


-3.48 


-11.88 


-20.11 


V 


0.82 


0.50 


0.39 


0.36 


0.35 


0.35 


0.34 


0.33 


0.31 


0.29 


s 


179.87 


168.27 


148.53 


139.70 


113.53 


99.67 


87.73 


77.52 


68.97 


62.10 


V 


0.99 


0.95 


0.81 


0.55 


0.25 


0.05 


0.21 


0.33 


0.43 


0.49 


6 


-1.31 


-4.79 


-6.79 


-6.48 


-4.03 


40.24 


67.02 


64.93 


60.66 


55.86 


V 


0.99 


0.93 


0.83 


0.72 


0.63 


0.54 


0.47 


0.41 


0.36 


0.31 


5 


1.63 


4.71 


4.03 


-0.24 


-6.81 


-14.75 


-23.54 


-32.88 


-42.63 


-52.76 


1] 


0.99 


0.93 


0.81 


0.68 


0.56 


0.48 


0.41 


0.35 


0.31 


0.27 


S 


-0.43 


-1.43 


-2.83 


-5.73 


-10.51 


-16.86 


-24.39 


-32.75 


-41.74 


-51.22 


V 


1.00 


0.99 


0.96 


0.89 


0.76 


0.59 


0.43 


0.34 


0.35 


0.40 


5 


-2.70 


-10.23 


-16.84 


-20.52 


-21.74 


-17.29 


-9.54 


2.48 


13.45 


19.05 


V 


0.98 


0.88 


0.73 


0.59 


0.49 


0.42 


0.36 


0.32 


0.28 


0.24 


S 


1.34 


5.24 


7.05 


5.13 


0.27 


-6.53 


-14.57 


-23.42 


-32.84 


-42.76 


V 


1.00 


1.00 


0.98 


0.93 


0.84 


0.73 


0.60 


0.48 


0.38 


0.31 


s 


0.13 


1.42 


3.34 


4.89 


5.42 


4.40 


1.51 


-3.28 


-9.93 


-18.44 


rj 


1.00 


1.00 


0.98 


0.94 


0.85 


0.73 


0.62 


0.52 


0.44 


0.37 


S 


-0.04 


-0.34 


-0.52 


-0.38 


-0.47 


-1.67 


-4.53 


-9.07 


-15.06 


-22.19 


1] 


1.00 


1.00 


1.00 


0.99 


0.98 


0.95 


0.91 


0.83 


0.73 


0.62 


5 


-0.20 


-2.28 


-5.76 


-9.30 


-12.17 


-13.97 


-14.52 


-13.70 


-11.56 


-8.11 


V 


1.00 


1.00 


0.99 


0.95 


0.88 


0.79 


0.69 


0.60 


0.52 


0.46 


s 


0.08 


0.99 


3.00 


5.76 


8.54 


10.44 


10.92 


9.91 


7.65 


4.44 


V 


1.00 


1.00 


1.00 


1.00 


0.99 


0.97 


0.92 


0.83 


0.71 


0.58 


s 


0.01 


0.31 


1.11 


2.12 


3.15 


4.21 


5.23 


5.95 


5.85 


4.32 


V 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.97 


0.94 


0.90 


0.85 


5 


0.00 


-0.10 


-0.30 


-0.41 


-0.24 


0.31 


1.25 


2.44 


3.67 


4.69 


V 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.98 


0.96 


6 


-0.01 


-0.47 


-1.85 


-3.76 


-5.78 


-7.67 


-9.28 


-10.49 


-11.25 


-11.50 


rj 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.97 


0.93 


0.89 


0.83 


5 


0.00 


0.17 


0.71 


1.62 


2.86 


4.41 


6.19 


8.00 


9.59 


10.71 


V 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.98 


0.97 


5 


0.00 


0.06 


0.39 


0.90 


1.47 


2.03 


2.59 


3.17 


3.81 


4.55 


V 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.99 


6 


0.00 


-0.02 


-0.12 


-0.26 


-0.36 


-0.36 


-0.20 


0.14 


0.68 


1.39 


V 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


S 


0.00 


-0.09 


-0.61 


-1.53 


-2.69 


-3.92 


-5.13 


-6.26 


-7.26 


-8.10 


V 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.98 


s 


0.00 


0.03 


0.21 


0.56 


1.05 


1.66 


2.42 


3.30 


4.33 


5.45 
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TABLE IX. Phase-shift and inelasticity parameters of the coupled partial waves with J = 1,2 
assuming isospin symmetry. 6, ej and ojj are given in degrees. 





100 


200 


300 


400 


500 


000 


700 


800 


900 


1000 




Vs 


0.63 


0.37 


0.22 


0.15 


0.14 


0.14 


0.16 


0.17 


0.19 


0.22 




Ss 


160.16 


141.90 


125.45 


111.64 


99.16 


86.70 


76.21 


69.19 


65.54 


64.21 






3.71 


9.61 


15.24 


20.32 


24.93 


29.08 


31.90 


33.49 


34.89 


37.29 




£1 


1.62 


6.90 


12.60 


17.51 


20.74 


21.16 


19.42 


16.88 


14.46 


12.61 




Vd 


1.00 


1.00 


0.99 


0.98 


0.96 


0.92 


0.88 


0.84 


0.79 


0.74 




Sd 


0.17 


2.01 


5.83 


11.08 


16.81 


21.13 


22.53 


21.15 


17.91 


13.54 




VS 


0.63 


0.43 


0.32 


0.25 


0.20 


0.17 


0.14 


0.12 


0.13 


0.18 




Ss 


173.41 


166.53 


159.57 


152.23 


144.07 


134.41 


121.85 


104.94 


87.16 


76.85 






-0.91 


-1.57 


-1.60 


-1.83 


-3.08 


-6.20 


-12.04 


-20.87 


-30.83 


-41.24 




£1 


-0.71 


-2.87 


-5.09 


-7.00 


-8.64 


-10.00 


-10.48 


-8.07 


-1.29 


4.78 




Vd 


1.00 


1.00 


1.00 


0.99 


0.97 


0.94 


0.89 


0.83 


0.76 


0.69 




Sd 


-0.03 


-0.49 


-1.42 


-2.43 


-3.20 


-3.60 


-3.79 


-5.09 


-10.94 


-21.33 




VP 


0.96 


0.74 


0.42 


0.21 


0.10 


0.07 


0.06 


0.06 


0.05 


0.04 


13p2 


Sp 


0.70 


3.82 


6.18 


7.70 


14.21 


27.68 


35.57 


37.95 


41.92 


51.99 






-0.76 


-2.19 


-2.36 


-1.09 


2.34 


8.38 


13.50 


17.31 


20.64 


22.71 




£2 


-0.18 


-2.11 


-5.30 


-8.55 


-11.14 


-11.37 


-9.80 


-7.77 


-3.92 


3.71 




Vf 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.97 


0.94 


0.89 


0.84 




Sf 


0.01 


0.23 


0.82 


1.44 


1.68 


1.49 


2.06 


3.35 


5.59 


9.60 




VP 


0.99 


0.91 


0.75 


0.57 


0.42 


0.32 


0.26 


0.21 


0.17 


0.15 


33p2 


5p 


0.38 


2.34 


5.48 


8.41 


10.48 


11.72 


12.32 


12.79 


14.08 


17.35 




0J2 


0.28 


0.85 


1.01 


0.65 


-0.11 


-1.18 


-2.49 


-3.96 


-5.49 


-6.91 




£2 


0.06 


0.67 


1.74 


2.89 


4.01 


5.19 


6.58 


8.31 


10.38 


12.49 


33p2 


VF 


1.00 


1.00 


1.00 


1.00 


0.99 


0.97 


0.93 


0.85 


0.73 


0.58 




Sf 


0.00 


-0.06 


-0.29 


-0.64 


-1.00 


-1.26 


-1.40 


-1.63 


-2.48 


-5.12 
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TABLE X. Phase-shift and inelasticity parameters of the coupled partial waves with J = 3,4 
assuming isospin symmetry. 6, ej and uj are given in degrees. 



7nai,(MoV/r) 


100 


200 


300 


400 


500 


GOO 


700 


800 


900 


1000 






1.00 


1.00 


0.97 


0.89 


0.74 


0.58 


0.47 


0.42 


0.41 


0.41 




5d 


0.02 


0.47 


2.52 


7.07 


13.52 


19.85 


24.11 


24.78 


21.37 


14.52 






-0.24 


-1.04 


-1.71 


-1.68 


-0.63 


1.37 


4.05 


7.18 


10.92 


15.92 




£3 


-0.01 


-0.41 


-1.59 


-3.29 


-5.03 


-6.37 


-7.10 


-7.34 


-7.44 


-7.86 




VG 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.99 


0.96 




h 


0.00 


0.03 


0.23 


0.57 


0.94 


1.26 


1.61 


2.16 


3.07 


4.40 




VD 


1.00 


1.00 


1.00 


0.99 


0.96 


0.91 


0.84 


0.76 


0.68 


0.61 




6d 


0.00 


0.08 


0.56 


1.97 


4.68 


8.61 


13.05 


16.96 


19.53 


20.49 






0.10 


0.42 


0.71 


0.79 


0.57 


0.01 


-0.85 


-1.96 


-3.23 


-4.57 




£3 


0.00 


0.14 


0.51 


1.05 


1.67 


2.33 


3.04 


3.86 


4.92 


6.38 




VG 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.99 


0.97 




h 


0.00 


-0.01 


-0.07 


-0.20 


-0.37 


-0.56 


-0.71 


-0.79 


-0.71 


-0.35 




Vf 


1.00 


1.00 


1.00 


1.00 


0.99 


0.96 


0.92 


0.87 


0.82 


0.79 




5f 


0.00 


0.04 


0.29 


0.96 


2.25 


4.16 


6.50 


8.79 


10.48 


11.18 




U)4 


-0.10 


-0.45 


-0.93 


-1.23 


-1.20 


-0.76 


0.15 


1.59 


3.63 


6.36 




£4 


0.00 


-0.08 


-0.52 


-1.26 


-2.19 


-3.19 


-4.18 


-5.10 


-5.91 


-6.64 


"iJ4 


Vh 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 




Sh 


0.00 


0.00 


0.06 


0.21 


0.43 


0.70 


0.99 


1.27 


1.57 


1.88 




rjF 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.99 


0.98 


0.95 


33^4 


5f 


0.00 


0.00 


-0.01 


0.05 


0.29 


0.80 


1.62 


2.69 


3.86 


4.86 




U14 


-0.93 


-2.63 


-3.95 


-4.78 


-5.30 


-5.65 


-5.91 


-6.14 


-6.36 


-6.57 




£4 


0.00 


0.03 


0.17 


0.41 


0.72 


1.05 


1.41 


1.77 


2.15 


2.56 


33/^4 


Vh 


1.00 


1.00 


1.00 


1.00 
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FIG. 15. (Color online) The Argand diagrams for the uncoupled S, P, and D waves, assum- 
ing isospin symmetry. The symbols on the lines denote the values of the antiproton laboratory 
momenta. 



VIII. SUMMARY 

In summary, motivated by renewed experimental interest in low-energy antiproton-proton 
scattering, we have presented a new energy- dependent PWA of all pp scattering data below 
925 MeV/c antiproton laboratory momentum. We have improved the model independence 
and quality of the PWA by using for the long-range interaction, next to the electromag- 
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FIG. 16. (Color online) The Argand diagrams for the coupled S-D, P-F, and D-G waves, assum- 
ing isospin symmetry. The symbols on the lines denote the values of the antiproton laboratory 
momenta. 



netic potential, the charge-conjugated one- and two-pion exchange potential derived from 
the effective chiral Lagrangian of QCD. We have updated the database and included the 
high-quality differential cross sections and analyzing powers for charge-exchange scattering 
pp — !■ nn that were measured in the last years of operation of LEAR. The final database 
contains 3749 scattering data, which are fitted with an excellent Xmin/^dat = 1.000 or 
Xmin/^df = 1.048. This implies that the long-range potential provides an excellent descrip- 
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tion of pp elastic and charge-exchange scattering, which we count as a success for chiral 
effective field theory. Further improvement of the PWA is certainly possible, but it will 
require additional high-quality experimental data. Below 400 MeV/c, there are hardly scat- 
tering data available. Spin observables will further constrain the PWA solution, provided 
they are precise enough. The results presented in this paper will serve as the starting point 
for more specific investigations of low-energy antiproton-proton scattering. 
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APPENDIX 



In this Appendix, we study in more detail the statistical quality of the final antiproton- 
proton database, by investigating the distribution of the contributions of the N^at = 3749 
individual data points to the total x^? Xtot 25|. In the PWA, this distribution is given by 



Pi 



analysis 



1=1 



(44) 



In Fig. [17] we plot this distribution as a histogram and compare it to the theoretical 
distribution for 1 degree of freedom, 

1 



,^-l/2^-t/2 



(45) 



In order to make this comparison quantitative, we give the moments, the central moments, 
and the corresponding errors for the distributions. For a distribution P{t), with t > 0, we 
define the moments fi'^ and the central moments fin by 



dtp{t) r , 

POD 

fln= / dtP{t)it-ft[r , 

Jo 



respectively. The errors on the moments are given by 
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FIG. 17. (Color online) Probability distribution functions versus x^- The tails, with the values 
> 2, are enlarged by a factor of 20. The histogram contains 3749 data points in bins with 
Ax^ = 0.05. 

and similarly for a^^. The lowest moments and their errors are given in Table IXII The 
agreement between the moments of -Pi(x^) and -Pi,anaiysis(x^) is reasonable, but not perfect. 

In fact, for two reasons Pi(x^) is not the best distribution to compare to. First, while the 
first moment of Pi(x^) is n[ = 1, that of Pi,anaiysis(x^) is yu'i = xL/^dat- Since (xL) = ^df, 
we should compare to a narrower distribution P(x^) = P^^Piif^^^X^) with /3 = A^df/^dat- 
Second, the data points with individual Xi > 9 were rejected, which affects the tail of the 
distribution and the higher moments. Therefore, it is better to compare Pi,anaiysis(x^) to 

\x'r'^'e~^'/'-'9i9-x'), (48) 

where 'y{s,z) = f^^^e^^dt is the lower incomplete gamma function and a is a constant 
chosen to satisfy (x^) = A^df/^dat! in our case, A^^f = 3578 and Nda,t = 3749, therefore we 
have cr = 0.989 and 7(i, ^) = 1.768. The Heaviside step function 6^(9 — x^) removes the 
tail with x^ > 9. Pi,o-,cut(x^) is also plotted in Fig. [17] and its lowest moments with errors 

49 



-Pl,<T,CUt(X^) 



^,19 



TABLE XI. Moments /i^ and central moments fin of the database of the PWA and of the two 
theoretical probability distribution functions. The errors are given for A'^jat = 3749, where the 
contributions of the normalization data are included. 







^l,(T,CUt(X^) 


-^1, analysis (X ) 




1.00 ± 0.02 


0.95 ± 0.02 


1.00 ± 0.02 


^2 


3.00 ± 0.16 


2.59 ± 0.11 


2.80 ± 0.12 




15.0 ± 1.6 


10.7 ± 0.7 


11.8 ± 0.8 




105 ± 23 


56 ± 5 


62 ± 6 


^2 


2.00 ± 0.12 


1.67 ± 0.08 


1.80 ± 0.08 




8.0 ± 1.3 


5.1 ± 0.5 


5.4 ± 0.5 




60 ± 18 


26.8 ± 3.1 


28.5 ± 3.1 



are given in Table |XT] as well. The agreement between the moments of -Pi,anaiysis(x^) ^^^^ 
-Pi,cr,cut(x^) is good, which implies that the distribution of the PWA is close to what is 
expected for statistical data. 
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